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NOMENCLATURE 

a  Channel  half -height 

A  Cross-sectional  area  of  the  tube 

B  Magnetic  induction,  B  =  (xeH 

C  Friction  factor  defined  in  equation  (29),  Chapter  $« 

D  Tube  diameter 

De  Equivalent  diameter,  I4.  x  cross-sectional  area/perimeter 

e  Electric  field  magnitude  factor,  Eq  =  euB0 

E  Electric  field  strength 

F  Body  force  vector 

FX,F„,FZ   x,  y,  z  components  of  body  force 

Fr,F/j,Fx   r,  ft,    x  components  of  body  force 

gx  x  component  of  the  acceleration  of  gravity 

H  Magnetic  field  intensity 

K  Constant  defined  in  equation  (29),  Chapter  5 

J  Electric  current  density 

M  Hartmann  number,  M  =  u.eHoa(o~e  ffp)  V2 

p  Fluid  pressure 

Pq  Pressure  at  channel  mouth 

p„  Pressure  at  far  upstream  where  fluid  velocity  is  zero 

A  p  Total  pressure  drop  in  the  entrance  region 

(p  -  Po) 
P      Dimensionless  pressure,  — — 

r 
R      Dimensionless  radial  coordinate,  — 


/>ua 

Raa 

Reynolds 

number, 

H 

ReD 

Reynolds 

number, 

V- 
i^us 

Res 

Reynolds 

number, 

V- 

Rer 

Reynolds 

number, 

P^o 

rQ  Tube  radius 

r,jif,x  Cylindrical  coordinates 

s  Channel  height,  s  =  2a 

t  Time 

u,v,w  x,  y,  z  components  of  the  velocity  field 

u  Average  velocity 

Uq     Center  line  velocity 

u 
U      Dimensionless  velocity  in  x-direction,  — 

u 
U|.     Velocity  at  fully  developed  condition 

u-"-     Deviation  velocity 

r  jVjjjV     r,  0,  x  components  of  the  velocity  field 

V  Velocity  vector 

avP 

V  Dimensionless  velocity  in  y-direction,  

V- 

x,y,z   Rectangular  coordinates 

HX    (x/a) 
X      Dimensionless  x-coordinate, 


X'     Dimensionless  coordinate, 


,Vu  Re£ 

jj.x  X 


/>s2u       k 
Entrance   length 


CHAPTER  1.   INTRODUCTION 

When  a  viscous  fluid  enters  a  conduit  (a  circular  tube  or 
a  flat  duct)  the  velocity  profile  of  the  flow  changes  as  the  flow 
progresses  from  the  entrance  to  the  interior  of  the  conduit  be- 
cause the  wall  of  the  conduit  tends  to  retard  the  flow.   Even- 
tually, the  velocity  profile  will  assume  a  form  which  is  invar- 
iant with  respect  to  the  direction  of  flow.   The  governing 
equations  for  the  flow  in  this  region  of  changing  velocity  pro- 
file, the  entrance  region,  are  the  continuity  equation  and  the 
equation  of  motion  given  below  (1,  2). 

9  P  -* 

+  div  (/°V)  =  0  (1) 

dt 

and 

0  DV  _>   _,    _, 

r  —  =  -grad  p  -  (V  •  t)  +  P.  (2) 

Dt 

If  the  flow  is  incompressible,  the  above  systems  of  equa- 
tions are  simplified.   When  temperature  variations  are  small, 
the  viscosity  may  be  taken  to  be  constant.   For  Newtonian  flow 
equations  (1)  and  (2)  then  become 

div  V  =  0  (3) 

and 

a™  2  •*   -» 

r  —  =  -grad  p  +  nV   V  +  P.  (k) 

Dt 
In  rectangular  coordinates,  equations  (3)  and  (lj.)  assume 
the  form 


0u        2v        Jw 

—  +  —  +  —  =  0  (5) 

3 x       ,?y       ,?z 

.  5u  ^u  #u  ^u  ^p 

P( —  +   u  —  +v  —  +   w  — )    = 

2t  ?x  gy  ?z  2x 

3  2u       <?2u       ^2u 

+   n( +  —-  +  —-)    +  Fx  (6a) 

dy?        2y2       ^ 

3V  ,?v  3V  gv  yp 

P  { —   +   u  —   +   v  —  +   w  — )    = 

dt  3x  Jy  ?z  dy 

3  2v       ^2v      ^2v 

+   n( r  +  —z  +  — r)    +  F  (6b) 

,?x^       ^       ^  J 

,?w  gu  3w  jw  jp 

f{ —  +  u  —  +v  —  +  w  - — )    =   -  — 

dt  3x  !)y  3z  3z 

0  2w        32w       32w 

+  lx( — ; +  ~~S  +  T}  +  Fz  (6c) 

3*2      3y2      3z2 

where  u,    v,    and  w   are   velocity   components   and  F_,    F    ,    and  Fz 

are   body   forces   in  the   direction  of  x,    y,    and   z,    respectively. 

In   cylindrical   coordinates   they   are 

3vT       vr        1     $vd        3vx 

+  —  +  -     — L.  +  =  o  (7) 

3v  r        r     dp  3  7. 

„   ^vr  3vr        v0    ^vr        v02  3vr  3V 

P( +vr  +  —  — —  +   vx  =  Fr 

it  dv  r     d</>  Jr  3x  Jv 

32vr        1  3vv        vr        1      32vr        2      3vd        32v 

+   |i{— —  + ___  +  __  —_  -  -r  —f-  +  — — )       (8a) 

3v2  V    VT        TZ        v2    9<f>  v23$  3xd 

.  3v0  2vo(        V0   2v0        v  v^  av0  1  ^p 

/3( +   vr  —!-  +  -£  — £  +  C.  +   Vx  _£)    =  p 

<?t  Jr  r       2j#  r  ^x  ^        r  P0 


32vri        1     w        vd        1     J2v0        2     ^vr        92v0 

+  n( — "  +  — -  -  -4  + + +  — -)    (8b> 

<?  r2        r  ^r  r2        r23  02  r2    3  0        ,?x2 

*vx  3vx        vd   <?vx  3v  «?p 

H +  vr  +  -£•  —  +  vx  )    =  Fx  -   — 

3t  dv         v     3p  3x  $x 

32vx        1  gvx        1      ;?2vx       ;?2v 

^r2  r  #r  r^  20^  ^ 

where  v  ,  vd,  and  vx  are  velocity  components,  and  Fr,  F7,  and 

Px  are  body  forces  in  the  direction  of  r,  0,  and  x,  respectively. 

For  two-dimensional  flow  under  the  absence  of  body  forces, 

equations  (5)  and  (6)  become 

(9) 


:iOa) 


(10b) 

in  a  rectangular  coordinate  system.   (When  the  body  force  ex- 

ists,  replace  p  by  P  =  p  —  -/!• ,  where  F  =  V-H. , ) 

When  the  flow  is  axisymmetrical  the  corresponding  equations 

in  cylindrical  coordinates  are 

3vr        vr        &vK. 

+  —  +  ——  =  0  (11) 

^r  r  £x 

#vr  svr  gvr  1  2P 

2t  r  3v  *  gx  Par 


<?u       2v 

—  +  —   =  0 

3x       3j 

3u            #u 
—  +  u  —  + 
it            dx 

V 

5U 

<?y 

lap           ,?2u     <?2u 
-  -  —  +   V  (— -  +  — 5) 

/?  2>x              «?x^       ^y 

2v            ?v 
—  +   u   + 

2>t           ^x 

V 

1   ?P              02v       ^2v 

+  v  ( —  +  - — ) 

C  37              d*2       JJ2 

J2Vr  1     ^V  Vr  22Y 


^r2  r  ^r  r2        <jx2 


5vx      ^vx      ?vx     1  jp 

+  vr  +vx  =  -  -  — 

?t       3r       jx  r  <?x 


/a2vx   1  ,?vx   ,?2v 

+  V  (— —  + +  — r)   .  (12b) 

gvd  r  gv         ax- 

It  is  often  considered  that  the  change  of  the  velocity  pro- 
file is  caused  by  the  formation  of  a  boundary  layer  whose  thick- 
ness increases  to  the  radius  of  the  tube  toward  the  downstream 
direction.   Then  the  following  Prandtl's  boundary  layer  assump- 
tion can  be  used  to  simplify  the  equation  of  motion. 

22u     i        S2vx\  <?2u  (      32vr 

11)   — ~  I  or  — —  /  is  very  small  compared  with  or  

3t?      V   <?x2  /  3J2    \       3r2 

so  it  can  be  neglected. 

(2)  The  transverse  velocity  v  (or  vp)  is  small  in  comparison 

with  u  ( or  v_)  . 

,  dv 

(3)  Consequently  the  pressure  gradient  —  is  a  function  of  x 

dx 
alone. 

Therefore  equations  (10)  and  (12)  are  reduced  to  the  fol- 
lowing boundary  layer  equations.   In  rectangular  coordinates 

du     ju.  3u     1  dp     3  \x 

—  +  u  —  +v  —  = +  V  (13) 

3t     3x  $y  f   dx     5y2 

and  in  cylindrical  coordinates 

■?vx      d^x  <5vx     1  dp     a2M  1  ^vx 

+  vr  +  vx  ..-__+  V  (——  + 1),       (Ik) 

dt  d?  <?x      P   dx     dv2  r  3r 

If  the  flow  is  steady,  the  terms  —  and  are  zero  and 

3t  at 

the  governing  equations  in  the  rectangular  coordinate  systems 


#u       <?v 

—  +  —  =   0 

ax     aj 

U    +    V    = 

2x            ?y 

1   dp            ,?2u 
?  dx            ^y2 

(15) 
(16) 


The  corresponding  equations  in  the  cylindrical  coordinate 
system  are 

^+Il+^=0  (17) 

Jr  r  ^x 

0vx  ,5>vx  1   dp  ^2vx        1   ?vx 

vr  —  +  vx  —  = +  y  (— — -  + )  .        (18) 

dr  dx  f   dx  3vl  r  3v 

A  study  of  the  velocity  field  and  friction  factor  in  the 
entrance  region  of  tubes  and  ducts  is  of  practical  importance 
and  has  been  a  subject  of  investigation  for  many  years. 

In  this  entrance  region  we  wish  to  obtain  expressions  for 
(1)  the  pressure  drop  between  any  two  sections,  (2)  the  velocity 
distribution  at  any  section,  and  (3)  the  value  of  x  at  which  the 
fully  developed  flow  is  attained. 

The  two  simultaneous  equations,  equations  (15)  and  (16)  or 
equations  (17)  and  (18)  with  appropriate  boundary  conditions, 
are  sufficient  to  solve  the  two' unknowns,  u  and  v  (or  vx  and 
v_) .   The  pressure  ceased  to  be  an  unknown  function  since  it 
can  now  be  evaluated  from  the  potential  flow  in  the  central  core 
by  the  aid  of  the  Bernoulli  equation. 

Although  equations  (15)  and  (16)  or  equations  (17)  and  (18) 
are  much  simpler  than  the  original  equations  representing  govern- 
ing equations,  their  exact  analytical  solutions  have  not  been 
found  so  far.   However,  various  approximation  methods  of 


solution  have  been  employed  to  solve  the  entrance  region  prob- 
lem.  The  approximate  solutions  may  be  classified  in  four  gen- 
eral categories—the  momentum  integral  method,  linearization 
method,  matching  method,  and  finite  difference  method. 

In  the  momentum  integral  method  the  flow  is  divided  into  a 
boundary  layer  part  near  the  wall  and  a  potential  flow  part  in 
the  central  core.   A  parabolic  velocity  profile  (or  any  other 
similar  velocity  profile)  is  assumed  in  the  boundary  layer  and 
is  joined  with  the  center  core  velocity  profile  which  is  assumed 
to  be  a  straight  line.   A  momentum  integral  equation  is  derived 
based  on  the  momentum  conservation  principle.   This  approach  was 
devised  and  applied  by  Schiller  (3)  for  flow  in  a  circular  tube 
and  is  similar  to  the  Karman-Pohlhausen  momentum  integral  method 
applied  to  a  flat  plate. 

In  the  second  category  of  solutions,  the  inertia  terms  of 
boundary  layer  equations  are  linearized.   This  category  of  solu- 
tions is  capable  of  providing  continuous  solutions  for  the 
velocity  of  distribution  and  pressure  drop  in  the  entrance  region. 
This  class  of  solutions  for  circular  tubes  is  mainly  due  to  the 
work  of  Langhaar  (I4.)  and  Targ  (5).   Langhaar's  approach  was  ex- 
tended by  Han  (6)  to  obtain  solutions  for  rectangular  ducts, 
and  Targ's  approach  was  used  by  Lundgren  and  Sparrow  et  al  (7, 
8)  to  obtain  solutions  for  channels  of  parallel  plates  and  for 
ducts  of  arbitrary  cross  section. 

The  third  group  of  solutions  is  constructed  by  matching 
together  the  boundary  layer  solutions  which  are  valid  near  the 
entrance  and  the  perturbations  of  the  fully  developed  solutions 


which  are  valid  for  downstream.   This  class  of  solutions  was 
originally  given  by  Schlichting  (1,  9)  for  a  parallel-plate 
channel  and  was  also  applied  to  circular  tubes  by  Atkinson- 
Goldstein  (10)  . 

The  fourth  approach  involves  reduction  of  the  continuity 
and  momentum  equations  to  finite  difference  equations  which  are 
solved  numerically  on  an  electronic  digital  computer.   This 
method  has  been  applied  by  Bodoia  and  Osterle  (11),  Hwang  and 
Pan  (12),  and  Wang  and  Longwell  (13)  for  flow  in  a  parallel- 
plate  channel.   However,  it  should  be  indicated  that  the  solu- 
tions in  References  (11,  12)  were  obtained  by  solving  the 
boundary  layer  equation,  equations  (l£)  and  (16),  while  the  re- 
sults in  Reference  (13)  were  obtained  from  the  solution  of  the 
continuity  equation  and  the  Navier-Stokes  equation,  equations 
(9),  (10a),  and  (10b),  in  which  the  time  dependent  terms  are 
omitted  for  the  steady-state  condition. 

The  purpose  of  this  report  is  to  present  in  detail  the 
solutions  of  the  various  methods  discussed  above.   Moreover,  the 
finite  difference  method  is  applied  to  solve  the  entrance  region 
problem  of  MHD  flow  in  a  parallel-plate  channel  with  a  non- 
magnetically  fully  developed  velocity  profile  at  the  entry. 
The  results  are  then  compared  with  those  obtained  by  means  of 
the  momentum  integral  method  (II4.). 
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CHAPTER  2.   SCHILLER'S  METHOD 

Schiller's  method  is  also  known  as  Karman-Pohlhausen 
momentum-integral  method.   When  a  viscous  fluid  enters  a  circu- 
lar tube  with  a  uniform  velocity  profile  across  the  cross  sec- 
tion, a  thin  boundary  layer  is  often  considered  to  have  been 
formed.   This  is  mentioned  in  the  previous  chapter.   The  thick- 
ness of  this  layer  increases  with  the  distance  'from  the  entry 
until  it  is  equal  to  the  radius  of  the  pipe.   Outside  the 
boundary  layer,  the  flow  is  considered  to  be  frictionless  and 
accelerated.   In  this  method  the  velocity  profile  is  assumed  to 
be  a  curve  build-up  by  two  parabolas  and  a  straight  line  (see 
Pig.  1) ;  the  vertices  of  the  parabolas  lie  on  the  border  of  the 
boundary  layer.   The  parabolic  velocity  in  the  boundary  layer 
is  assumed  to  be 

u     y     7  ? 

—  =  2(-)  -  (-r  CD 

u0      6      5 

where  y  =  rQ  -  r 

and  Uq  is  the  velocity  of  the  central  core  and  is  a  function  of 

x;  6  is  the  thickness  of  the  boundary  layer. 

The  equation  of  continuity  and  the  equation  of  motion  (the 
momentum  equation)  give  two  relations  among  u0,  x,  and  6.   Elim- 
ination of  6  will  lead  to  a  relation  between  u0  and  x. 

The  momentum  theory  is  used  here  for  the  x  direction  to 
obtain  a  momentum  integral  equation  representing  the  equation 
of  motion.   The  control  surface  is  formed  by  two  planes  a-b  and 
c-d,  spaced  apart  from  each  other  by  the  small  distance  dx  and 
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by  the  wall  (see  Pig.  2).   By  means  of  the  momentum  balance  on 
this  control  surface  the  following  momentum  integral  equation 
is  obtained  (see  Appendix  I  for  detail). 

—  [        u2(r0  -  y)dy   -  u0  —  J  u(r0  -  y) dy 
dx  Vq  /       dx  V0 


in     J-    „  /   o      u   p\    qua 
-2Vuo(^)  +-ro2f2_-  (-)2)u0-^ 


ro 


dx 


12) 


Substituting  the  assumed  velocity  profile  expressed  by 
equation  (1)  into  the  integral  form  of  continuity  equation  and 
carrying  out  the  integration,  we  obtain  a  relation  between  the 
boundary  layer  thickness,  5,  and  the  center  line  velocity,  uQ , 
as  shown  below: 


5 

rn 


(3) 


J+  -  6(1 ) 

L'o  uo 

where  u  is  the  average  velocity. 

Substituting  equations  (1)  and  (3)  into  equation  (2),  and 

carrying  out  the  integration,  we  obtain  a  relation  between  uQ 

and  x  as  follows: 

x/D       1  [58     22  17   ; 

=  —  j  —  y\ in(l  +  l\) yi).  +  2W  -  2Vlz 

16  I  1^    5  15 


Re 


D 


16,  A  -  2)7 

5 


37 


;'2 


i  +  n 


10 


.   2H  -  1 

I1  (— — ) 


37l/2~ 


L0 
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Therefore  the  velocity  profile  in  the  entrance  region  is  defined 
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by  equations  (1),  (3),  and  (1|)  .   Given  a  value  of  uQ,  the  values 
of  x  and  6  are  known  from  equations  (Ij.)  and  (3).   The  velocity- 
distribution  at  the  corresponding  section  is  defined  by  equation 
(1).   The  flow  is  fully  developed  when  6  =  rQ  and  u0  =  2u. 
Hence  Y\    -   1.   Then  from  equation  (L\.)  ,    the  entrance  length  is 
found  to  be 

xe  =  0.0288  DReD. 

Differentiating  both  sides  of  equation  (I4.)  with  respect  to 

d  >j           du0 
x,  we  see  that  ,  that  is  ,  does  not  vanish  at  Y\    =1. 

dx  dx 

This  means  that,  the  center  line  velocity  does  not  asymptotically 

approach  it3  final  value,  2u,  at  the  fully  developed  section. 

For  the  calculation  of  pressure  drop  in  the  entrance 

region,  the  Bernoulli  equation  is  applied  since  it  is  assumed 

that  a  potential  flow  in  the  central  core  exists  throughout  the 

entrance  region.   Therefore  the  pressure  distribution  is  assumed 

to  be  dependent  only  on  the  potential  flow.   The  pressure  drop 

between  any  two  locations  in  this  region  is  found  as 

Pi  -  P2  =  ~  (<u0)22  "  ("O)!2]  =  -  ^   (  (122  "  V) 

+  2(n2  -  vj  • 

Since  x  =  0  and  '1  0  =  0  at  the  entry,  and  x  =  xe  and  11  =  1  at 
the  fully  developed  section,  the  total  pressure  drop  in  the 
entrance  region  is 

Ap  =  -  Pu2    . 
2 

The  pressure  drop  in  the  entrance  region  is  larger  than 
that  in  the  fully  developed  region  of  equal  length  because  of 


15 


the  increased  friction  loss  end  the  fact  that  the  kinetic  energy 
of  the  fluid  increases  as  it  passes  downstream.   For  the  evalua- 
tion of  total  pressure  drop  between  a  section  in  the  entrance 
region  and  a  section  in  the  fully  developed  region,  it  is  more 
convenient  to  treat  the  entrance  region  as  a  fully  developed  one 
and  the  pressure  drop  thus  calculated  is  modified  by  a  correc- 
tion factor.   The  equation  for  calculating  the  total  pressure 
drop  in  the  fully  developed  section  is  found  to  be 
p0  -  p    61+ 


P u2    ReDD 


x  +  1.16 


2 

where  1.16  is  the  correction  factor.   If  the  fluid  enters  the 

pipe  from  a  region  where  the  pressure  is  x>00     and  the  velocity 

is  negligible,  there  is  a  pressure  drop  of  —  f\i,    and  then  the 

2 
total  pressure  drop  is 

p   -  p   6kx 

=  +  2.16  . 


f>\x2  ReDD 


2 
The  derivation  and  solution  are  given  in  detail  in  Appendix  I. 
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CHAPTER  3-   LANGHAAR 'S  LINEARIZING 
APPROXIMATE  METHOD 


The  difficulty  in  obtaining  an  exact  analytical  solution 
of  the  boundary  layer  equations  stems  from  the  fact  that  the 
convective  terms  are  nonlinear.   Langhaar  (1)  introduced  an 
assumption  that  the  convective  terms  in  the  boundary  layer  equa- 
tion for  the  axial  direction,  x,  (see  Fig.  1)  are  equal  to  the 
product  of  the  kinematic  viscosity  \>,    the  axial  velocity  u,  and 
a  parameter  [3.   This  parameter  is  considered  to  be  a  function 
of  x  only.   Langhaar 's  assumption  can  be  written  as 

?U     5U      ?U       2 

u  —  +  v  —  +  w  —  =V(3u  (1) 

?x     jy     3z 

and  gives  rise  to  the  following  linearized  form  of  the  boundary 

layer  equation 


32,. 


32n       g   u 


Vp^u  =  gx  -  -  —  +  J>  (—  +  — )  (2) 

|°  dx      3j2       d^ 

where  gx  is  the  x  component  of  the  acceleration  of  gravity. 

In  order  to  solve  the  problem  for  a  flow  in  a  circular 

tube,  equation  (2)  is  transformed  into  an  expression  with 

cylindrical  coordinates  as 

d2u   1  du 

__  + p2u  =  ^  (3) 

dr^   r  dr 
where 

1       1  dp 

V  f>    dX 

which  is  a  function  of  x  alone. 
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Equation  (3)  is  regarded  as  a  velocity  function  generating 
equation.   The  solution  of  equation  (3)  is  expressed  in  terms 
of  a  modified  Bessel  function  as 

i0(/)  -  lot/—) 

u  r0 

-  = (U 

u       I2(/) 

where  7*  =  (3r0,  rQ  is  the  radius  of  the  tube,  and  u  is  the 
average  velocity. 

Therefore  if  the  relation  between  f    and  x  is  known,  equa- 
tion (I4.)  will  define  the  velocity  profiles  everywhere  in  the 
entrance  region.   It  is  worth  noting  that  when  (3  =  0,  (that  is, 
7—0),    the  acceleration  terms  are  zero.   The  solution  then 
assumes  the  form  of  a  fully  developed  parabolic  profile.   When 
P  — ^  co  ,     (that  is,  / — 5.  co)  ,    the  solution  assumes  a  uniform  pro- 
file at  the  entry  (see  Appendix  II  for  details).   Thus  the  so- 
lution provides  a  smooth  transition  from  the  entry  where 
y  — ■>  <*£>    to  the  fully  developed  section  where  /  =  0. 

For  determining  function  "/(x),  the  momentum  integral  equa- 


tion is  used.   Integrating  the  boundary  layer  equation  in  cylin- 

v/e 

drical  coordinates  over  the  cross  section,Aobtain: 

a  fr°   9  _   1  dp   „      (r0   ,2U   i  au 

2  —    u^rdr  =  gxr</ r^  +  2  V  j   (— -  + )rdr.  (5) 

dx  >Q  P   dx         /Q   3rd        r  3v 

2   1  dp 
The  term  gxrQ   -  —  —  rQ  can  be  eliminated  by  applying  the 
r   dx 

Navier-Stokes  equation  in  the  center  line  of  the  tube.   After 
the  substitution  of  equation  (I4.)  and  simplification,  the  follow- 
ing relation  is  obtained. 
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A/2 

X  =      YdZ  (6) 

h 


where 


x 

X  =  — 


r0   Rer 

I2(/) 
Y  =  -y 


and 


27ix{f)  -r2 

l4.I0(r)l2(/)  -  [lo(/)  -  l]2  -  2Ii2(/) 


2I22(/) 
Equation  (6)  provides  a  relation  between  x  and  Y.      The  inte- 
gration can  be  carried  out  numerically;  the  results  are  pre- 
sented in  Table  1  of  Reference  (1). 

The  entrance  length  is  defined  as  the  distance  from  the 
entry  to  a  section  where  the  center  line  velocity  is  99  per  cent 
of  the  fully  developed  one.   This  definition  differs  from  that 
given  by  Schiller.   When  the  value  of  99  per  cent  is  reached, 
the  ratio  of  the  center  line  velocity  to  the  average  velocity  is 
1.98.   From  equation  (1+)  , 
Uq  _  Iq(Y)    -  1 

u      I2(Y) 

u0 
Hence,  when  —  =  1.98,  the  corresponding  value  of  /(and  there- 

u 
fore  the  value  of  x)  can  be  determined.   The  entrance  length 

is  found  to  be 

•xe  =  0.0568  D  ReD  . 

In  'Calculating  the  pressure  drop  in  the  entrance  region, 

Langhaar  applied  the  mechanical  energy  balance  instead  of  using 

the  Bernoulli  equation  because  in  the  latter  part  of  the  entrance 
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length,  the  Bernoulli  equation  ceases  to  be  applicable.   The 
energy  loss  due  to  viscous  dissipation  was  considered.   Equat- 
ing the  work  of  pressure  force  and  frictional  force  to  the 
change  in  kinetic  energy  gives  the  energy  equation 

/X     I         0U.    o        1    f  _r> 

uA(p0  -  p)  =  (i,  I  dx  /   ( — )"   dA  +  -f>   \     (u3  -  UJ)  dA.    (7) 
'0    'A        dr  2   lA 

Substituting  equation  (I4.)  into  equation  (7)  and  simplifying, 

gives  the  total  pressure  drop  between  a  supply  reservoir  and 

the  section  x; 

X  r2  A 

— 2  (I12  -  I0I2)dX  +  2  j   (U3  -  l)RdR   (8) 
0   x2  '0 


p<=<,  -  p 

=   1   +   2 

2 

where 

u 

U    =  - 

u 

and   R  =  —  . 
r0 
The  integrations  in  equation  (8)  can  be  carried  out  numerically 

and  the  results  are  presented  in  Table  2  of  Reference  (1) . 
Hence  the  pressure  drops  in  the  entrance  region  are  defined. 

To  calculate  the  pressure  drop  between  the  entrance  of  the 
tube  and  a  location  in  the  fully  developed  region,  the  entrance 
region  is  treated  as  if  it  were  a  fully  developed  one  and  the 
pressure  drop  thus  calculated  is  modified  by  a  correction 
factor  as  follows: 
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+  2.2£ 


1    „    DRen 
-P   u2  D 

2 

where  2.28  is  the  correction  factor. 

The  solution  in  detail  is  given  in  Appendix  II. 

Reference 


H.  L.  Langhaar,  "Steady  Flow  in  the  Transition  Length  of 
a  Straight  Tube,"  Trans,  of  ASME,  J.  of  Applied  Mechanics, 
6k,    A -55,  (191^2);  also  see  Ph.D.  Thesis  by  H.  L.  Langhaar, 
Lehigh  University,  (1914-0). 
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CHAPTER  l+.       SCHLICHTING'S  MATCHING-  METHOD 

Let  us  consider  the  flow  in  a  parallel  plate  channel  with 
a  uniform  velocity  distribution  at  the  entry  (Reference  1) . 
It  is  assumed  that  the  flow  is  incompressible  and  two-dimen- 
sional.  The  governing  boundary  layer  equation  and  continuity 
equation  are 

Ju     ,?u     1  dp     32u 

u  —  +  v  —  = +  V  — ;  (D 

3u   2v 

—  +  —  =  0  .  (2) 

In  this  method  the  flow  field  is  divided  into  two  sections, 
the  upstream  and  the  downstream.   In  the  upstream  section,  the 
boundary  layers  develop  in  a  similar  way  to  the  flow  over  a  flat 
plate  with  a  pressure  gradient  in  the  direction  of  flow.   The 
stream  function  and  center  line  velocity  are  expanded  in  power 
series,  as  in  the  case  of  the  flat  plate.   The  solution  is  found 
after  the  coefficients  of  the  series  are  determined  by  substi- 
tuting the  assumed  series  into  equation  (1).   For  the  downstream 
section,  the  velocity  is  assumed  to  be  the  sum  of  the  fully  de- 
veloped parabolic  velocity  distribution  and  a  deviation  velocity 
which  approaches  zero  as  x — ?  oo.      This  deviation  velocity  is 
obtained  by  a  series  expansion  in  the  upstream  direction.   Hav- 
ing obtained  both  solutions  in  the  form  of  a  series  expansion, 
they  are  joined  at  a  point  where  both  solutions  are  valid.   In 
this  way  an  approximate  description  of  the  flow  field  in  the 
entire  entrance  region  is  obtained. 
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For  the  solution  of  the  upstream  section,  the  center  line 
velocity  is  assumed  to  be  in  the  following  form 

u0(x)  =  u(l  +  Kx  c  +  K2  e2  +  .  .  .  )  (3) 


»x 
where   £  =J        _    ;    and  K1;  Ko  .  .  .  are  constants  to  be  de- 

"  a^  u 
termined. 

Substituting  equation  (3)  into  the  pressure  gradient  at 
the  center  line  in  the  flow  direction,  we  obtain 
1  dp      du0   u2£  _ 
f   dx      dx     2x 

+  3(K1K2  +  K3)£   +  .  .  .J     (Ij.) 
It  is  also  assumed  that  the  stream  function,  f(x,  y'),  can  be 
expanded  into  a  series  in  £  .   This  series  is  similar  in  form 
to  that  of  u0  but  with  coefficients  which  are  functions  of  x  and  y' 
as  shown  below. 

-2, 


V(x,   y')   =  ua[ef0C>l)    +  e2f1(7j)   +   .    .    .] 


where  1  =  h/ J ' 

y  x 


Then 


2  f  - 

u(x,    y<)    = =  u[f0'(?|)    +     £fx<(Tj)    +    .    .    .  "jf         (5) 


v  ( x,    y ' )    =   - 


- —  =  —  r^f0  +  2  £2fi  +  3  £2f2  + . . .] 

3X  2x  l  ^  J 


ua  £A| 

+  


[f0'    +     6^1+     62f2>    +   .    .    .].    (6) 


2x 

Inserting  equations  (I).),  (5),  and  (6)  into  equation  (1)  and  col- 
lecting terms  in  powers  of  6  ,  we  obtain  differential  equations 
for  f0,  f1,  ...   .   The  equation  of  the  first  approximation 


23 


(for  the  zero  power  of  c  )  is 

f0f0"  +  2  f0"'  =  0  .  (7) 

This  is  identical  to  the  differential  equation  of  flow  over  a 
flat  plate  given  and  solved  by  Blasius.   The  second  approxima- 
tion (for  the  first  power  of  <=■  )  is  determined  by  the  equation 

2f1'"  +  fofx"  -  f0'fi'  +  2  rQ"fl   =  -K-l  (8) 

which  is  a  linear  differential  equation  of  third  order.   The 
boundary  conditions  are  f^  =  f-,  '  =  0  at  7\  =  0   and  f -,  '  =  "K-. 
at  7]  =oO, 

These  differential  equations  can  be  solved  numerically. 
With  the  values  of  f„,  f,,  .  .  .,  and  K-,,  K2,  .  .  .  being 
found,  the  solution  of  this  upstream  section  is  complete. 

For  the  solution  of  the  downstream  section,  we  assume 

u(x,  y)  =  uf(y)  -  u-::-(x,  y)  (9) 

where 

3_    J2 
Uf(y)  =  -  u(i -) 

2        a2 
and  u-»-  is  the  deviation  velocity,  which  is  assumed  to  be  in 
the  form 

—  =  cl6  1        f  1    (-)  .  (10) 

u  a 

Substituting  equations  (9)  and  (10)  into  equation  (1)  yields 

the  following  differential  equation  for  <f ,  '. 

rf  +  3^g  (i  -^>vl+  vj  =  o  (id 

with  X-^  as  its  eigenvalue.   The  boundary  conditions  are 


2k 


^  =  V 

=   0 

at 

y  =  0 

fl  =  ?l' 

=   0 

at 

y  =  ±a 

Vx'  -  i 

at 

y  =  o   . 

The  last  boundary  condition  is  chosen  arbitrarily.   This  is 
possible  since  C]_  is  still  free.   The  solution  for  <f,  in  equa- 
tion (11)  is  obtained  by  the  power  series  method.   Assume 


y   S?    y 


^(-J  =  L     An(-)n  .  (12) 

a    n=0    a 


By  substituting  equations  (12)  into  (11),  the  recursion  formula 
for  An  is  obtained. 

_  3X1[(n-2)(n+l)An  -  (n-1)  (n+2) An+2] 
U+k  2(n+l)(n+2)(n+3)(n+W 

The  values  of  An  and  \^  can  be  obtained  from  this  recursion 
formula  together  with  boundary  conditions.   Thus  the  deviation 
velocity  is  defined  if  Cx   is  known,  which  can  be  determined  in 
the  process  of  matching  both  solutions. 

The  pressure  drop  in  the  upstream  section  is  obtained  from 
Euler's  equation 

1  dp      dug 

"  3  —  =  u0  •  (13) 

I     dx      dx 

Since  the  boundary  layer  thickness  is  very  thin  in  this  section 
and  the  flow  is  mostly  irrotational  in  the  central  core.   Inte- 
grating equation  (13)  gives  the  pressure  drop  between  any  two 
locations  in  the  upstream  section.   On  the  other  hand,  the 
pressure  drop  in  the  downstream  section  involves  the  perturbation 
velocity,  u-::-(x,  y)  .   Substituting  the  assumed  velocity  distribu- 
tion, equation  (9),  into  the  boundary  layer  equation  (l)and 
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noting  that  the  pressure  varies  with  x  only,  we  obtain 

VV"  (l)J  .  (Ik) 


i  f£ = v  r  —  -  —  e"xie2 

r  dx     I-  a  2    a2 


Integrating  equation  ( 11+)  from  the  joining  section,  £  ■::-  to  £ 
gives  the  pressure  drop  'in  the  downstream  section 


2 

For  the  complete  determination  of  the  velocity  distribu- 
tions as  well  as  the  pressure  drops  in  the  entire  entrance  region, 
only  C-,  and  £*,  the  location  of  the  matching  section,  remain  to 
be  determined.   This  is  done  when  the  upstream  and  downstream 
solutions  are  matched.   The  conditions  that  the  center  line 

velocities  of  both  solutions  be  equal  and  that  the  pressure 

dp 
gradient  —  be  identical  at  the  matching  section  provide  two 

dx 
simultaneous  equations  for  C-,    and  6*  as  follows: 

3  2 

1  +  K-l  6  +  K2£2  +  I^e3  +  K^  =  -  -  Cie"  1  ,  (15) 

(1  +  K-^  +  K2&2  +  K3£3  +  Kr£^)  (I^  +  2Kg£  +  3^2   +ljKi]£3) 

=  3  +  Cl9"Al   fy"(l)  .       (16) 
Therefore  C]_  and  £*•  are  determined  from  equations  (15)  and  (16). 
The  solution  is  then  complete. 

The  numerical  values  of  C-l  and  ~k-^   given  by  Schlichting  are 
0.3t|.85  and  18.75,  respectively.   Values  of  Kx,  K2,  ...,  V1,   f1", 
...,  and  Tq,    f]_,  ...,  are  not  available.   The  values  calculated 
by  others  can  be  found  in  References  (2,  3) • 
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The  entrance  length  given  by  Schliohting  Is 


xe  =  0.010  HeD  De 


where  De  is  the  equivalent  diameter. 

The  solution  is  given  in  detail  in  Appendix  III. 


References 

1-    YorkCh}960ting'iAQ,UndarX  LaJer   The0r?'  McGraw-Hill,  New 

lork   I960,  p   169,  or  H.  Schlichting,  "Laminare  Kanalein 
laufstromung,  "  ZAMM  11;,  pp.  368,  (1934)  -         ^analein- 

2.    M.  Roidt and  R.  D.  Cess,  "An  Approximate  Analysis  of 

ofT^latTc't^f  °„df~TS°-P^  *?  thS  EntraLe^^ion 
PP.  171-176  Applied  Mechanics,  March,  1962, 

3'    mJl'll'+fFl    "MaSnetohy<arodynamic  Effect  upon  Laminar 
p  ™  "*  Entrance  Region  of  Perfectly  Conducting 
Parallel  Plates,"  M.S.  Thesis,  University  of  Pittsburgh, 


27 


CHAPTER  5.   A  FINITE  DIFFERENCE  ANALYSIS  OF 

MAGNETO HYDRODYNAMICS  CHANNEL  ENTRANCE 

FLOW  WITH  PARABOLIC  VELOCITY 

AT  THE  ENTRY 


Summa  ry 

Tne  entrance  region,  problem  of  magnetohydrodynamics  flow 
with  nonmagnetically  fully  developed  (parabolic)  velocity  at 
the  entry  is  investigated.   The  flow  is  in  a  horizontal  channel 
of  parallel  and  electrically  nonconducting  plates  on  which  a 
uniform  magnetic  field  is  imposed  vertically.   The  basic  equa- 
tions are  transformed  to  finite  difference  equations  and  solved 
by  a  digital  computer.   The  variations  of  velocity  profiles, 
friction  factors,  and  pressure  drops  for  Hartmann  numbers  of 
2.5,  L\.,    10,  and  $0    are  presented.   Values  of  the  entrance 
length  are  compared  with  the  only  available  results  by 
Maciulaitis  and  Loeffler,  who  solved  the  problem  by  momentum- 
integral  method.   The  two  sets  of  results  appear  to  be  sub- 
stantially different. 

Introduction 

This  chapter  presents  a  study  of  magnetohydrodynamic  (MHD) 
flow  in  the  entrance  region  of  a  parallel-plate  channel.   The 
velocity  profile  at  the  entry  of  a  MHD  channel  is  assumed  to  be 
a  fully  developed  non-MHD  flow  (parabolic  profile) .   The  velocity 
profile  eventually  develops  to  be  a  Hartmann  velocity  profile 


28 


in  the  entrance  region  of  a  parallel-plate  channel  where  a  uni- 
formly distributed  magnetic  field  is  applied  to  the  channel. 
The  entrance  length  of  such  a  MHD  channel,  the  development  of 
the  velocity  profile,  and  the  pressure  drop  are  studied  in  this 
work. 

Although  the  analysis  of  MHD  entrance  region  problem  has 
been  made  by  many  investigators,  most  of  them  (1,  2,  3>  h-) 
studied  the  case  of  developing  velocity  from  a  uniformly  dis- 
tributed one  at  the  entry  to  Hartmann  velocity.   There  is  only 
one  paper  (5)  published  studying  the  case  of  development  of 
velocity  from  parabolic  to  Hartmann  velocity.   In  this  paper  the 
Karman-Pohlhausen  integral  method  with  the  assumption  of  Schil- 
ler's velocity  profile  was  used.   Since  there  is  no  other  liter- 
ature available  dealing  with  this  specific  case,  the  results 
could  not  be  compared,  although  the  results  obtained  for  the 
case  of  uniform  to  Hartmann  velocity  by  using  the  same  integral 
method  have  agreed  well  with  others  (5) ■ 

A  practical  importance  in  this  problem  is  the  length  of  the 
entrance  region.   The  length  of  the  entrance  region  is  defined 
in  Reference  (5)  as  the  distance  required  for  the  friction 
factor  to  come  within  10  per  cent  of  the  fully  developed  value. 
Another  usual  definition  of  the  length  of  entrance  region  is 
the  length  required  for  the  center  line  velocity  to  reach  99 
per  cent  of  the  fully  developed  one.   The  results  obtained  in 
the  present  work  are  quite  different  from  each  other  using  the 
different  definitions.   Also  the  result  in  the  present  work  is 
quite  different  from  that  of  Reference  (5)  using  the  friction 
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factor  definition.   The  entrance  length,  defined  according  to  the 
comparison  with  center  line  velocity  is  not  available  in  Refer- 
ence (5)  • 

A  finite  difference  analysis  method  (2,  3)  which  yields 
good  results  for  the  developing  velocity  from  uniform  to  Hart- 
mann  velocity  is  used  in  the  present  work. 

Basic  Equation 

Let  us  consider- an  electrically  conducting  fluid  entering 
a  horizontal  channel  of  parallel  plates  which  are  electrically 
nonconducting.   The  velocity  profile  has  been  assumed  to  be 
nonmagnetically  fully  developed  (parabolic  profile)  at  the  entry 
of  the  duct  (that  is,  x  =  0),  after  which  a  uniformly  distributed 
magnetic  field  is  applied  perpendicularly  to  the  plates.   (See 
Fig.  3-)   The  following  assumptions  are  made  for  simplicity: 
(a)  the  flow  is  laminar  and  steady;  (b)  all  fluid  properties 
(f ,    Cp,  k,   ])  )    are  constant;  (c)  permeability,  p.e,    and  electrical 
conductivity,  5~e,    are  constant  scalar  quantities;  (d)  rapid 
oscillations  do  not  exist, and  therefore  the  displacement  current 
is  negligible;  (e)  the  effect  of  gravitational  force  is  negli- 
gible; (f)  the  usual  Prandtl  boundary  layer  assumptions  apply; 
(g)  variation  in  the  z-direction  is  negligible,  i.e.,  the  flow 
is  two-dimensional;  (h)  the  electric  field  term  measured  across 
the  insulated  channel  walls  is  zero;  and  (i)  the  induced  field 
in  x-direction,  Bx,  is  negligible  in  comparison  with  the  applied 
field,  Bq , 
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The  governing  basic  equations  may  be  written  as  follows  (1) 
Maxwell's  equations  in  mks  units  are 


(1) 
(2) 

(3) 


curl  H  =  J,    .,  div  J  =  o 

curl  E  =  -ix  —  ,        div  H  =  0 

Ohm's  law  for  moving  fluids  is 

J  =  o-     (E  +  V  x  u  H) 
-*  e 

where  V  is  total  velocity  in  vector  form. 

The  continuity  equation  is 

div  V  =  0  . 

The  modified  Navier-Stokes  equation  is 

3~^  1 

-+  (V.  graa^  =  --gradp+l/v2^  +  _  {J  ,  ^)  .         {S) 

Under  the  assumptions  made  above,  the  basic  equations  (1) 
to  (5)  can  be  simplified  to  the  following  two  equations: 

9U  .    3U     1  dp     d2u  B 

In  the  potential  flow,  equation  (7)  may  be  written 
^   duo  _    1  dp    <reB02       <reB0E0 
°  dx    "  ?  dx  "    />   U°      T~ 
1  dp    tf"eBQEo  _     du0    <TeBo2 

"?dx"    ^>   =U°^r  +  -7T-uo  (8) 

where  u0  is  the  center  line  velocity. 

Substituting  equation  (8)  into  (7),  we  have 


or 


-    V^V  +  T(U°-U)+P37  •   (9) 

It  may  be  need  that  the  electric  field,  E0,  does  not  appear  in 
equation  (9)  and  thus  has  no  effect  on  the  velocity  profile. 
However,  from  equation  (8),  it  does  have  an  influence  on  the 
pressure  gradient  required  to  maintain  the  velocity  distribution 
in  equation  (9) . 

The  maximum  value  for  EQ  is  obtained  by  assuming  that  the 
side  plates  are  open-circuited.   This  permits  maximum  build-up 
of  the  electric  field  and  is  equivalent  to  no  net  current  in 
the  z-direction,  or 


a 


Jzdy  =  o  .  (10) 

Substituting  Jz  =  re(S0   +  uB0)  into  equation  (10),  one  obtains 

J       <re(E0  +  uB0)dy  =  creE02a  +  <reB0  /   udy  =  0  .    (n) 

'-a 
Since  the  continuity  equation  can  be  written  as 


/•"' 


j       ^  =  2Ua  (12) 

'-a 

substituting  equation  (12)  into  (11)  results  in 

H0(max)  =  ~u  B0  . 
Introducing  the  electric  field  factor,  e,  which  varies  betwe 
zero  and  one  with  the  external  resistance  varying  from  zero  t 
infinity,  E0  may  be  written  as 

E°  =  '™   Bo  •  (13) 

Combining  equations  (13)  and  (7)  gives 


en 

o 
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2u  2u  i   dp  02u        cre302 

u  _  +   v  _  =--  —  +JJ    — -  +  — — -   (eu   -  u)     .         (llj.) 
5x  ^y  Z5  dx  Py2  /> 

Introducing   the   following  ditnensionless   parameters 

i-ix  x/a 


?a2u 

Rea 

Y  =   y/a 

U  =  u/u 

av/1 
V   — 

1* 

P  _  p  "  po 

/>u2 

7~' 

M  =   tieK0a^ 

/  (Te 

7v  ' 

i    Hartmann  number, 

equations    (6),     (12),    and    (ll|.)    become 
3V        $V 


+   —   =   0 

ax     n 

r1 

TJdY  =   1 

'0 

2U             2?U 

U   —   +   V  —  = 

PX            3Y 

dp        ,?2U 
dX       <?Y2 

(15) 
(16) 

+  M2(e    -  U)     .  (17) 


The  boundary  conditions  on  the  flow  field  are: 
at     X  =  0,  for     -1  <■  Y   <  1  : 

U  =  1.5(1  -  Y2),    V  =  0,    P  =  P0  =  0 

for    X  >  0  and     Y  =  0: 

ra 

—  =0,      V  =  0 
dY 


(18) 
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for    X  >  0  and     X  =  +1: 

U  =  0,     V  =  0  . 

The  first  boundary  condition  shows  that  the  velocity  profile  is 
parabolic  at  the  entry  (X  =  0)  after  which  a  uniformly  dis- 
tributed magnetic  field  is  applied. 

Solution  of  the  Equation 


The  method  applied  follows  closely  that  of  Hwang  and 
Fan  (3)  .   Considering  the  mesh  network  of  Fig.  I4.,  a  grid  is  im- 
posed on  the  flow  field.   Equation  (17)  is  then  written  in 
implicit  finite  difference  form  at  any  point  (j,  k)  as  follows: 


U 


j,k 


u]+l,k  '  uj,k 


AX 


J,k 


uj,k+l  -  uj,k-l 
2AY 


3+1 


P  A  1    U 


j+l,k+l  -  2uj+l,k  +  Uj+l,k-] 


AX 


2  L 


uj,k+l  "  2uj,k  +  Uj^k-l" 


(AY)2 

9,    uj+l,k  +  Uj.k 
+  M2(e  -  — — ) 


(AY)2        J  (19) 

Considering  the  U  and  P  values  in  the  j+1  column  as  unknowns, 
and  those  in  the  j  column  as  known  values,  an  equation  for  U  and 
P  is  obtained. 

<-<k>  ^j+i,k-i  +  (Pk)  uj+i)k  +  (yk)  uj+1)k+1  +  (|)  pj+1  =  (j«k) 

(20) 


wnere 

■*k  = 

Pk  = 

1 

1 

_2(£  Y)2_ 

r u^ 

Ax 

M2 
2 

(AY)2 

n/L- 


n-l 


\ 


fe+1 


3p 


'a  v : 


•.  +  i 


T 


5      jv. 


j,!:+:     ; -;•;..;*.; 


i         I         I 


i*I.h 


o         i 
=  o 


-V' f 


-i         j  jfl 


Flq.4.     Mesh    r.atwork    for  -divvc-ronce 


representation. 
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£     = 


r 


.2(AY)2  J 


1 
~Ax 


J 


j,k 


AX 
1 
2(AY)2 


■     V-!     V      t J       "     — 

J'  2AY  AX 

(Uj,k-1    -    2Uj)k   +Uj>k_!)    -   M2(< 


UJ,< 


Using  Simpson's   integration  rule   for   equation    (l6),  tve  obiain 

1  11-0 

UdY   =  -    .    (Ux   +   ^U2    +    2U3   +   IjITjl    +    2U£   +    .     .     .)    =   1 


Ul   +   ^U2k   +   2  ^    U2k+1   =   3n' 


l,    2,   3, 


(21) 


Now,  considering  the  points  k  =  1,  2,  .  .  . ,  n  in  the  j+1  column 

in  equation  (20),  and  combining  them  with  equation  (21),  we 

obtain  (n  +  1)  equations  and  the  same  number  of  unknowns.   Here 

wo  use  the  boundary  conditions  that  —  =0  for  x  2  0  and  Y  -  0 

3Y 


and  that  U 


j+l,n+l 


yields  Uj+i,2  -  Uj+1,0 
1,  2,  .  .  .   . 

The  equations  in  matrix  form  are  as  follows. 


=  0  for  j  +  1  =  0, 
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1  k 

2 

J4- 

0!   ■<!+'/! 

0 

0 

-<2        ?2 

Y2 

0 

o        -^ 

?3 

''3 

2 1+  0 

o o  y 

0 0  % 

0 0  | 


■*n-l  Pn-1   ?n-l   £ 


uj+l,l 

r  3n  -1 

UJ+l/2 

0i 

uj+l,3 

02 

Via 

?3 

Vl,H 

Ci 

Vl        J 

A. 

The  values  of  U,  P,  and  V  are  known  at  the  section  X  =  0  from 
boundary  conditions.   Thus  let  j  =  0  at  this  section;  then  the 
values  of  U  and  P  in  the  j  +  1  =  1  column  can  be  obtained  by- 
solving  the  (n  +  1)  linear  simultaneous  equations.   In  order  to 
solve  V  in  this  column,  equation  (15)  is  replaced  by  the  follow- 
ing difference  equation: 

Pj+l,k+l  +  uj+l,k  -  Pj,k+1  -  uj,k   vj+l,k+l  -  vj+l,k 


2ZiX 


Ay 


Solving  for  Vj+1  k+1,  one  obtains 

AY 
Vj+l,k+l  =  Vj+1)k  -  —  (Uj+i)k+1  +  Uj+1>k  -  uj;k+1  -  Uj;k) 

(23) 
Thus  let  k  =  1,  2,  .,.,    n  and  substitute  the  values  of  U  which 
have  been  calculated  into  equation  (23),  and  remember  V -+1  ,  =  0 
from  boundary  conditions;  then  values  of  \f-j_  k  (k  =  1,  2,  ...,  n) 
are  obtained.   After  the  values  of  U,  V,  and  P  are  determined 
for  the  j  -1=1  column,  the  above  procedure  is  applied  to 
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the  j  +  1  =  2  column.   Thus  it  is  possible  to  advance  column  by 
column  along  the  channel  until  the  flow  becomes  fully  developed. 
The  (n  +  1)  simultaneous  equations  were  solved  by  Gaussian 
elimination  method;  the  calculations  were  performed  on  an  IBM 
1620  digital  computer.   For  Hartmann  numbers  equal  to   2.5,  k, 
and  10,  a  1^0-point  mesh  was  used  with  AX  =  0.00025  from  X  =  0 
to  X  =  0.001.   Then  a  20-point  mesh  with  AX  =  0.001  was  used 
until  fully  developed.   Because  of  the  relatively  shorter 
entrance  length  and  in  an  attempt  to  obtain  better  results  for 
the  case  where  Hartmann  number  equals  $0 ,    an  80-point  mesh  with 

Ax  =  0.000,05  is  used  from  X  =  0  to  X  =  0.0001.   Then  a  I4.O- 
point  mesh  with  AX  =  0.0003  is  used  until  X  =  0.001.   Thereafter 

AX  is  changed  to  0.00025  with  l4.0-po.lnt  mesh  until  fully  de- 
veloped.  Owing  to  the  capacity  of  the  digital  computer,  two 

programs  are  used  successively  in  the  case  of  80-point  mesh. 

AX  1  1 

In  all  cases,  the  ratio  x  was  kept  ■<  -  ,  and  as  close  to  - 

(AY)2  2  2 

as  possible  (3).   Each  of  the  above  listed  meshes  (80,  lj.0,  20 

points)  applied  to  one-half  the  channel  height. 

Results  and  Discussion 

The  velocity  component,  U,  for  Hartmann  number,  M,  of  If, 
10,  and  50  in  the  entrance  region  is  shown  in  Figs.  5a,  5b,  and 
5c  for  various  values  of  X.   It  is  the  same  as  in  the  case  of 
developing  velocity  from  a  uniform  velocity  distribution  at  the 
entry,  the'  dimensionless  center  line  velocity  approaches  1.281). 
for  M  =  l+,    1.111  for  M  =  10,  and  1.020  for  M  =  50 .   These  values 
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■  =  r. 


1+2 


can  be  obtained  from  the  velocity  equation  of  Hartmann  (6), 
which  is 

cosh  M  -  cosh  (MY) 


U  =  M 


[2k) 


_M  cosh  M  -  sinh  M  . 
The  velocity  profile  for  a  different  Hartmann  number,  M,  in  the 
entrance  region  is  shown  in  Pigs.  6a  and  6b.   It  is  well  known 
that  the  effect  of  a  magnetic  field  is  to  flatten  the  velocity 
profile.   Therefore  flows  with  velocity  profiles  nonmagnetically 
fully  developed  (parabolic  profile)  at  the  entry,  in  general 
require  larger  entrance  length  than  flows  with  uniform  velocity 
profile  at  the  entry.   Table  1  shows  the  entrance  length  for 
one  per  cent  and  five  per  cent  deviation  of  center  line  velocity 
from  fully  developed  one,  and  are  compared  with  values  from 
Hwang  and  Pan's  (3)  analysis  of  flows  with  a  uniform  entry 
profile. 

Maciulaitis  and  Loeffler  (5)  defined  the  entrance  length  as 
the  distance  required  for  the  friction  factor  to  come  within  10 
per  cent  of  the  final  developed  value.   In  order  to  make  com- 
parisons with  their  results,  the  entrance  length  according  to 
this  definition  for  different  M  wea  also  found.   By  definition, 
the  friction  factor  f  is 

_  tw    _  -Ou/;?y)y=a 


Pu2/2  fu2/2 

In  dimensionless  parameters,  it  becomes 


2   3U 
Rea  n 


Y=l  ' 


(25) 


(26) 


Therefore. the  ratio  of  the  friction  factor  to  that  at  the  fully 
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Table  1.   Comparison  of  the  entrance  length  based 
on  center  line  velocity. 


K  =  1+ 


M  =  10 


Xe'(l#) 
Xe'(l#) 

xe>(#) 


:  Parabolic  entry  :  Uniform  entry 
:  (Present  work)   :   (Ref.  (3)) 


0.019,3 

0.008, 18 

0.007,31 
0.003,77 


0.018,8 
0.007,32 

0.003,04 
0.000,234 


M  =  5o 


Xe'(l#) 
Xe'(^) 


0.000,275 

0.000,162 


developed  condition  is 


LP.D. 


<?U 
(— ) 
?Y  Y=l 


9U 


3V 
(—) 
3Y   Y=l  at  P.D. 


(27) 


The  value  of  ( — )     at  the  fully  developed  condition  is  ob- 

3Y  Y=l 
tained  from  equation  (2l\.)    as  follows. 


?U 

(— ) 
<?Y  Y=l  at  F.D. 


Mc  sinh  M 


(28) 


M  cosh  M  -  sinh  M 
Table  2  shows  the  comparison  of  the  present  work  with  that  of 
Maciulaitis  and  Loeffler  employing  the  entrance  length  in  accord- 
ance with  their  definition  for  M  =  If.   and  10.   Because  of  the 
drastic  change  of  velocity  near  the  wall  for  M  =  50,  the  fric- 
tion factor  obtained  by  the  finite  difference  technique  is  not 
accurate  due  to  the  capacity  of  the  available  computer.   There- 
fore only  M  =  I4.  and  10  are  shown  in  this  table. 


M> 


Table  2.   Comparison  of  the  entrance  length  based 
on  friction  factor. 


d„„„„„4-  ,.„„i,  i-r    \\         '■      Maciulaitis  & 
Present  work  (X  '  )      r— «n«.«.  tv    i\ 


e  '    :   Loeffler  (Xa' ) 


M  =  2.5  0.0038 

M  =  k  0.0045  0.0223 

M  =  10  O.OOI67  0.00^96 

The  disagreement  between  the  present  work  and  Maciulaitis 
and  Loeffler1 s  is  seen  in  Table  2.   In  their  solutions,  there 
is  a  minimum  value  of  6/a  for  certain  M,  for  example,  for  M  =  if, 
(6/a)min  is  0.612  (=1/0/11),  and  for  M  =  10,  (o/a)min  =  0 .  2l|495 
(=y6/l0).   Beyond  these  minimum  values  of  6/a,  the  solutions 
are  not  applicable.   Thus  the  velocity  profiles  corresponding 
to  these  (6/a)  j   become  the  fully  developed  one  which  was  used 
to  calculate  the  friction  factor  ratio.   It  is  doubtful  that  the 
fully  developed  velocity  profiles  thus  obtained  approximate  the 
exact  Hartmann  velocity  profile.   However,  the  Hartmann  velocity 
profiles  are  used  in  the  present  work,  and  the  velocity  dis- 
tributions in  the  later  part  of  the  entrance  region  are  very 
close  to  the  fully  developed  Hartmann  velocity  profiles  as  can 
be  seen  from  Pig.  6b. 

The  comparison  of  the  entrance  lengths  in  Tables  1  and  2 
shows  that  quite  a  large  difference  exists  using  different  defi- 
nitions.  It  is  debatable  that  the  definition  by  Reference  (5) 
can  be  used. 

The  accuracy  of  present  analysis  can  be  seen  in  a  reasonable 


k-1 


agreement  of  the  results  with  the  known  asymptotes,  for  instance, 

SV 
the  center-  line  velocities  and  —  evaluated  at  the  walls,  as 

3Y 
shown  in  Table  3. 


Table  3.   Comparison  of  the  asymptotic  values. 


Asymptotic 
Center  line  velocity 


Expected: Pre  sent: Deviation 
:work   :    {%) 


3U 
Asymptotic  ( — ) 

3X   Y=l 


Expected: Pre sent: Deviation 
:work   :    (%) 


k  1.282+2  1.2862  0.156 
10  1.1110  1.1123  0.117 
50   1.0201+   1.0218    0.137 


5-3286   5-2563   1.356 
11.1111   10.3810   6.571 


In  general,  as  the  Hartmann  number  increases  the  entrance 
length  becomes  shorter,  as  shown  in  Table  1.   However,  for  very 
small  Hartmann  numbers,  as  the  Hartmann  number  decreases  the 
entrance  length  decreases  also.   In  particular,  when  Hartmann 
number  is  zero,  since  there  is  no  change  in  velocity  profile 
throughout  the  flow  field,  the  entrance  length  is  zero.   The 
same  conclusion  was  also  given  in  Reference  (5).   Table  1+  shows 
the  variation  of  entrance  length  defined  by  one  per  cent  devia- 
tion of  center  line  velocity,  with  M  of  1.5,  2.5,  1+,  10,  and  50 . 
The  maximum  entrance  length  will  occur  for  the  Hartmann  number 
of  approximately  2.5/~'lj-. 

The  dimensionless  pressure  drops  over  the  entrance  region 

1 

are  shown  in  Pig.  7  for  M  =  1+,  10,  and  50  for  e  =  -  .   Figure  7 

2 

is  not  convenient  for  calculating  the  total  pressure  drop  beyond 


^8 


z/ft  I 


rJ7 


k9 


Table  /j.. 


1.5 


2.5 


10 


5o 


Xe'(l$)     0.0145   0.0211   0.0205   0.0076   0.000,2£ 


the  entrance  length.   For  this  purpose,  the  following  equation 
is  used  to  express  the  pressure  drop  for  X1  greater  than  the 
entrance  length. 
Po  "  P 


Pu2/2 


=  CX'  +  K 


(29) 


where  C  and  K  are  constants  to  be  determined. 

Considering  the  fully  developed  region,  equation  (17) 
becomes 


dP   d2U 

—  =  +  M2(e  -  U)   . 

dX   dY2 

From  equation  (21j_),  we  have 

dU  -sinh  MY 

—  =  M2( )  , 

dY      M  cosh  M  -  sinh  M 

d2U  -cosh  MY 

=  M3{ )  . 

dY^      M  cosh  M  -  sinh  M 

Substituting  equations  (21;)  and  (31)  into  (30)  yields 


dP 


—  m2 


M  cosh  M 


M*(i 


-) 


dx  M  cosh  M  -  sinh  M 

Integrating  equation  (32)  from  Xe  to  X,  one  obtains 

M  cosh  M 


p  M  cosh  M       fX 

dP  =  M2(e )     dX 

p  M  cosh  M  -  sinh  M  /jj 

Afi  e 


? 


(30) 


(3D 


(32) 


5o 


M  cosh  M 
P  .  px  =  M2(e )  (X  .  t  ) 

M  cosh  M  -  sinh  M 

*»  -  ",  M  cosh  M 

^  =  2M2(e  -  )  (X  -  XJ 


u2/2  M  cosh  M  -  sinh  M 

The  total  pressure  drop  is  the  sum  of  the  pressure  drop  in  the 
entrance  region  and  the  pressure  drop  in  the  region  after  the 
velocity  is  fully  developed.   Therefore  the  total  pressure 
drop  is 

PQ    -   P  _  PO    '   PX9        Pxe    "   P 


pu2/2  f>u2/2  f>u2/2 

=   -2M2(e 


M    cosh  M  Pn    -    Py 


M   cosh  M   -    sinh  M 


pn    -   p  0  M   cosh  M 

-   -2M2(e   - 


/"uya  M   cosh  M   -    sinh  m' 

M  cosh  M 


+  2M2(e 


•"•e' 

r 

l2/2      ' 

X  + 

PO 

-  Px9 
:2/2 

)  xe 

(33) 

M  cosh  M  -  sinh  M 
Since  X  =  l\.X<  ,    equation  (33)  becomes 

p0  -  P      p         M  cosh  M  PO  "  Py 

±^-7-  =   -8M2(e )  X'  +  -   ,  Xe 

fu-^/2  M  cosh  M  -  sinh  M       /3u2/2 

M  cosh  M 

+  2M2(e )  X„  .         (3k) 

M  cosh  M  -  sinh  M 

Comparing  equation  (3I4.)  with  (29),  we  have 

M  cosh  M 

C  =  -8M2(e )  (35) 

M  cosh  M  -  sinh  M 
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K  = 


Po  -  PXF 


2M2(( 


M  cosh  M 


-)  Xfi 


(36) 


P\i2/2  M  cosh  M  -  sinh  M 

Table  5  shows  the  comparison  of  values  of  C  calculated  from 
equation  (35)  with  those  obtained  from  finite  difference  calcu- 
lation, and  also  shows  values  of  K  obtained  from  this  analysis. 
Their  excellent  agreement  provides  an  additional  measure  of 
confidence  in  the  results  of  the  present  analysis. 


Table  5. 


c 

M 

:      Prom   equation 

^5) 

Present 
work 

De 

viation 

:        K 

k 

106.628 

106.596 

0.03 

-0.232 

10 

I4.88.889 

tj.89.9ti4 

0.216 

-0.369 

50 

10l|0  8.162 

lOlj.36.896 

0.276 

-0.538 

It  is  obvious  from  equation  (36)  that  the  correction 

factor  K  is  a  function  of  e.   The  values  shown  in  Table  5  for 

1 
e  =  —  ,  and  Xe  is  the  entrance  length  listed  in  Table  1  (one 

2 
per  cent  deviation) .   It  is  interesting  to  note  that  the  values 

of  K  are  negative.   This  means  that  contrary  to  the  case  of 

uniform  velocity  profile  at  the  entry,  the  pressure  drop  in  the 

entrance  region  is  smaller  than  the  pressure  drop  in  the  fully 

developed  region  of  equal  length.   This  is  true  because  the  free 

stream  velocity  is  decreasing  as  it  passes  downstream.   This 

causes  the  pressure  to  increase  and  compensates  the  pressure 

drop  due  to  friction  in  the  entrance  region.   Thus  the  overall 

pressure  drop  becomes  smaller. 
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APPEHDIX  I.   SCHILLER'S  METHOD 

The  flow  Is  in  a  horizontal  circular  tube  with  radius  rQ . 
In  Schiller's  method  (1,  2)  the  velocity  profile  is  approxi- 
mated by  a  curve  comprised  of  two  parabolas  and  a  straight  line 
(see  Pig.  1);  the  vertices  of  the  parabolas  lie  on  the  border 
of  the  boundary  layer.   The  parabolic  velocity  profile  in  the 
boundary  layer  is  assumed  to  be 

_  =  2  (-)  -  (-)*  (1) 

Uq       6      6 

where  Uq  is  the  velocity  of  the  center  core  and  is  a  function 
of  x,  the  coordinate  along  the  direction  of  the  flow;  6  is  the 
thickness  of  the  boundary  layer,  and  y  =  rQ  -  r. 

The  flow  is  considered  to  be  steady,  laminar,  incompres- 
sible, and  axi symmetrical.   Furthermore,  it  is  assumed  that  the 

fluid  properties  are  constant,  that  the  pressure  is  uniform  at 

2P 
any  given  section  and  is  a  function  of  x  only  ( —  =  0),  and 

37 
that  the  velocity  at  the  inlet  edge  is  uniform. 

First  of  all,  the  momentum  integral  equation  for  flow  in  a 

circular  tube  is  derived.   Consider  the  flow  region  within  the 

boundary  layer  between  sections  ab  and  cd,  which  are  dx  apart. 

(See  Fig.  2.)   Let  the  x-axis  be  in  the  same  direction  as  the 

flow,  and  let  the  y-axis  be  measured  from  the  wall.   Then  the 

mass  flow  rate  into  the  section  ab  within  the  boundary  layer  is 

/°u(2itr)dr  . 


Since  r  =  r0  -  y,  y 
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r,  and  dy  =  -dr,  the  above  equation 


becomes 


f° 

Min  =   j     2*MrQ  -  y)(-dy) 


=  2it/>  I   u(rQ  -  y)dy  . 
'0 
The  mass  flow  rate  out  of  the  region  is 

/5  d 

Mout  =  27lH    u<ro  -  y)dy  +  2*P  — 

'o  dx 

The  difference  between  M_,,+  and  Mj_  is 

oat      m 


u(r0  -  y)dy 


L'O 


dx 


"out  "  Min  =  **f 


ax 


u(rQ  -  y)dy 


dx 


which  is  the  mass  flow  rate  into  the  region  from  the  center  core. 
Therefore  the  momentum  flux  into  the  region  due  to  the  bulk 
flow  (or  convective  flow)  is 


2%P  u0  — 
dx 


j  u(r0  -  y)dy 


dx  +  2nf       u2(rQ  -  y)dy 
0 


The  momentum  flux  out  of  the  region  due  to  the  bulk  flow  is 

,6 


^6  a 

2%f>  j       u2(r0    -   y)dy  +   271/3  — 


'0  d*  L'o 

Thus  the  net  momentum  flux  out  of  the  region  is 


u  (r0  "  y)dy 


dx 


r/6 


2%f 


dx  j_/0 


L0 


y)dy 


rva 


dx  -  2t[^  uq 


dx 


u(rQ  -  y)dy|dx.  (2 


For  steady-state  flow,  the  net  momentum  flux  out  of  the  region 
is  equal  to  the  net  force  acting  on  the  region  in  the  same  di- 
rection.  The  forces  exerted  on  the  region  are  the  pressure 
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force  and  the  shear  force  of  the  wall.   Hence 


Ipx  =   -Tw(2itr0clx)    +   p|*r02    -   Tt(r0    -   6)2] 

dp    r   ? 

+  (p  +  £  —  dx)  htr1   -  n(r1  -  46) 


dp 
(p  +  —  dx) 
dx 


nr02  -  n(r0  -  6  -  d6)2 


where  0  <  £  <  1.   The  first  term  of  the  right-hand  side  repre- 
sents the  shear  force  of  the  wall,  the  second  term  the  force 
acting  on  the  boundary  layer  at  section  ab,  the  third  term  the 
force  acting  on  the  boundary  from  the  central  core,  and  the  last 
term  the  force  acting  on  the  boundary  layer  at  section  cd. 
Rearranging  the  equation  yields 
■Z?x   =  -Tw(2itr0dx)  +  pn(2r06  -  62) 

+  (p  +  t   —  dx)  X  |2r,de  -  (d5)2 
dx       L 

-    (p   +  —  dx)    it|(2r06    -   62)    +   2rQd6    -    26d6    -    (d6)2| 
dx 

Neglecting  the   infinitesimal   terms   of  higher   orders  yields 

£Fx  =   -tw(2nr0dx)    +   Pt(2r1d6)    -   pn(2rQd6    -    26d5) 

a     dp 
-it(2r06    -   8*)    —   dx 
dx 


or 


ZFX  =   -Tw(2nr0dx)    +   pn(2r1d5)    -   p7i[2(r0    -   6)  do] 

o      dP 
-    7i(2r06    -    5^)    —   dx    .  (3) 

dx 


Since 

rn    -   6   =   r 


X, 


58 


pit[2(r0    -   6)ds]   =   pn(2r1d6)     . 
Prom   equations    (3)    and    (k) ,    wa    obtain 


(W 


dp 


Z.FX  =   -TwC2nr0dx)    -  n(2r06    -   6^)    —  dx 

dx 

Prom  Euler's  equation,  we  have 

dp       dug 

-  —  =  ho  —    ■ 

dx       dx 
From  the  definition  we  also  h8ve  (note  that  tw  is  the  stress 


exerted  by  the  fluid  on  the  wall) : 

'yo. 


fry 


Therefore 


^u                                    du0 
rPx  =   -2nr0   (i(— )      0   dx  +  npu0  (2rQ6    -   e^)dx   . 


dx 


Prom  equation    (1),    we    obtain 
3u  ua 

Substituting   the   above   equation   into   the   net   force   equation, 
we   obtain 

^                                    u0                                         o          du0 
HFX   =   -2ltr0    u-(2   — )dx  +   7t/>(2r05    -   6^)u0   dx 

6  dx 


r0  p 

=  -Ij.Ttnuo-C — )<3x  +  V'O 
5 


'6       6   o 
2(— )  -  (  — ) 
•  r0     r0   . 


dur 


dx 


According  to  the  momentum  tbeorom,  equation  (2)  is  set 
equal  to  the  right-hand  side  of  equation  (5)  bo  give 
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2%P- 


dx 


u^tr-Q    -   y)dy   dx  -    2-k/>uq   _ 
'0  dx 


1       u(rQ    -   y)dyldx 
fa  J 


-iwun( — )dx  +  jt/V. 
6  u 


-)    -    (- 


dun 


dx 


dx 


d 
dx 


u2(r0    -   y)dy 


=    -2H(-)    +_rQ2 


dx 


u(rQ    -   y)dy 


L'O 


5  6 

2 (— )2 

r0  r0 


dur 


a  x 


C6) 


Equation    (6)    is   a  momentum   integral   equation   of   the   boundary 
layer  for  flow   in   a   circular   tube.      Using   the   assumed   velocity 
profile   expressed   in  equation    (1),    we   can  perform  the   integra- 
tion  in   the   first   and    second   terms    on   the   left-hand   side   of 
equation    (6),    as  follows: 

r5  '     f5  y2        y3     y^ 

u2(r0    -   y)dy  =  u02(t)_  _  -  1*.  -_  +  __)     (r0    -   y)  dy 


=  u0  r0 


53        6^ 

'4    y3      it    y^      i    y5 


6^  3  6J  U.         6^  5 


k   ^     k   y5 


i    y 


6  -.6 


62   1|  63    5  6^   6     J0 


8  2        11     22 


u2(rQ    -   y)dy  =  u02(l6 

30 


6 
r0 


11  ) 

2 
r0 


(7) 


6o 


/5  y     y2 

u(r0  -  y)dy  =    /      u0(2 -)(r0-y)dy 

>0  5         5 


a0J0 


2   y2        1      y3l6 


6    2 
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2   y^        1      y^5 
6    3  62   Ij.   J0 


un5< 


-  ro5uo  -  —  u0 


2  a  K2 

r0  5  5 

u(r0   -  y)dy  =  u0(8  —   -  5  — - 

12  r0  r0 


(8) 


Substituting   equations    (7)    and    (8)    into   equation   (6),    we   obtain 


_2_  _ 
30      dx 


u02(l6 


11—) 

r0 


ro> 


=    -2"Vu0(— )    + 


r0<  d 

u0  — 

12  dx 

6  5? 


6  62 

V8-  -  5—) 

r0  vQd 
du0 


dx 


30    L 


u02  —   (16 11  )    +   2(16 11 


:)uf 


dx 
,2 


r0" 
6 


dx    J 


ro^  r        d  6  5' 

—  u0     u0  —   (8 5 = 

12  L        dx  r0  r02 


-)    +    (< 


-   V    (2   — 


r0" 

du0 


$■         dun 
"   5  — )   — 
r0        dx 


2  r0         r0 

Rearranging   the    terms,    we   obtain 

2 


o  -  -=  -2»V- 

dx  I 


,r0. 


r02             d                6  62  5  5J 
iV   —    (32 22   — -   -    ijfl   —  +   2^  — 


60 


dx 


r0 
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r„2    dun     5       62       6       62 
+  -2-  Uq  —2  (6^.  —  -  1+4  — -  -  ^0  —  +  25  — - 
oO      dx       r^        rvX        rn        -n^2 


r0       r0 


«*>' 


6o  —  +  30  — ) 


-2Vu0(— ) 

5 


—  uo  —  (3 
60     dx 


-)  + 


2  j       K2 

r0   duQ 


Tq     60   dx 


6 
(11  —  -  36  — ) 


=  -2vA 

6 


Simplifying  the  equation  yields 


u0  -  (3  — 
dx    rQ2 


L0 


duQ 
•)  +  —  (11 

dx 


62 


36 


120   r0 

—  y(-).  (9) 


<0 


r0 


If  u  is  the  average  velocity,  the  principle  of  conserva- 
tion of  mass  gives 
,r0 


jirQ2  u  =  (2nr)udr  +  u0it(r0    -   6)' 


rn-5 


=   2ti     I        (r0    -   y)    u(-dy)    +  uQ"K(r0    -   6)' 
'8 


=    2it 


J        (r0    -   y)    udy+u0T[(r0    -    5)' 


f  y      y 

2nu0  I     (r0    -   y)(2 -)  dy  +  u07t(r0    -   6): 

0  5        5 


2nr0u0 


2   y2        1      y3 


6    2  5^    3   J0 


r  —        "    2%un 


2  y3      i    yine 


6   3  5*=  I4. 


+  u0it(r0    -   5)' 
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or 

?   -         ^  ^               ?                                          9 

icrQ^  U  =  -  huqPqS  -  -  TtuQ6      +  im0(r0    -   5)   , 

3  6 


(  — )    -    -   (  — )2  +    (1 )2 


Therefore 

u  1+5  5     5 

u0        3      r0  6      r0  r0 

or 

U  2      5  1      5      0 

-  =   1    --(_)+   -    (_)2( 
u0  3      r0  6      r0 

Rearranging   the    equation,    we    obtain 

1   5   _   2   5  u 

-  (— )2  --(—)  +  (1 )  =  0 

6   ro     3   r0         u0 

or 

5   „      5           u 
(_)^  .  ^(  —  j  +6(1 )  =  0  . 

r0        r0  u0 

5 
This  is  a  quadratic  equation  of  ( — ) .   The  roots  are 

r0 

6        J  « 

-  -  2  +f  1+  -  6(1 )  . 

But 
5 

r0 

Hence 


Si 


5         A  u 

—  =  2  -  7.11.  -  6(1 )  .  do) 

r0  ^0 

Equation-  (10)  shows  the  relation  between  the  boundary  layer 
thickness  and  the  velocity  at  the  central  core. 
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jet    r. 


I   - 
u 


-  1 


then 


(11) 


UQ    1  +  7] 
Substituting  equation  (11)  into  equation  (10),  we  obtain 
6 


=  2  -1/1).  -  6(1 


1+7^ 


r0 


1+  -  2>1 


(12) 


■  1  +  1\ 

Prom  equation  (11),  we  can  obtain 

u0  =  (l  +  r|)u 
and 

du0   _  d  >| 

dx      dx 
Substituting  equations  (12),  (13),  and  (llj.)  into  equation  (9) 
results  in 


[13) 


(11).) 


(1  +  V])u  —  3(2 
dx 


It.  "  21 


1  +  1 


-)   -  8(2 


1  +  >1 


-) 


d  m 

+  u  

dx 

lj.80 


11(2  -  /)/ )   -  36(2 

1  +  r) 


1  +  V] 


2  - 


A  -  2  C] 


(15) 


1  +  n 

where   D  =  2  r- ,  the  diameter  of  the  tube. 

The  left-hand  side  of  this  equation  is  simplified  as 
follows. 


61* 


_   d  k~2n         4"2l  k.-2H 

u+nu  —  3(4  -  4v + )  -i6  +  8  y  — 

dx    L  1+1        '    1+  )\  1+  >] 


_  dn  /4-2i      4-2  n  /it-- 21 

+  u  11(4  -  4  V  +  )    -  72  +  36  1/  

dx  i+  n        i+  y\  i+ii 


d        12-6)1 

=    (l+rj)u  —   ( -   4 

dx        1+  V] 


4-21  _  Ah       44-22J] 

4)    +  u  ( 

1+  r\  dx  i+n 


4-2  n 


28) 


1+  r\ 
(l+^)(-2)-(4-2il 


(l+r])u 


(l+?])(-6)-(i2-6>7  )  l 

5 4  •  - 

(1  +  A)  2 


[i  +  n 


4  -  2n 
i  +n 


11 

dx 


_      d>,       44-22>j 

+  u    (— L)( 

dx  1+1 


d1 


dx 


-  d  1 

u  — - 

dx 


-18 


12 


1+7J  /4-21 

(i+n)- 


4-2^1 
i+  n 

44-22^1 
1   +  Y\ 


-    2t 


4-2  ?! 
1+1 


28 


-2-50^ 

.  i  +  n 


1+  >7 

12 


(1+11) 


4-2  1 

i+  n 


4-2  n 
i+  ^ 


Thus   equation    (15>)    becomes 

d^j    r-2   -   50  yj  12 

u  +  : 

dx     j       1  +  A 


(1+ *1) 


/Hi 

i+  a 


4-2  ^ 
i+1 


480 


2  - 


4-2  >i 
1+1 


65 


('2 


1+T\ 


-2  -  50 n 


12 


1  + 1 


l4_8o  V 

D      Du 


(i+»l! 


1+1 


14.-2*1 


11 

dx 


-1).   -    100^ 


2J+ 


-   16 


l  +  r\ 


$0T[ 


(i+y\) 


k-2n 

i+n 


k-zn 

1+  Y\ 


i  +  y\ 


12 


1+1 


1*1 


I+-27J 


1+  >] 


1  +  1 


l|8o     i 


D     Rer 


-1161  +   16 


21+ 


i  +  n 


5or? 


(i+n ) 


k-2  n 
i+ 1 


16 


k-2n 
i+  >i 


14--2  /? 
1+  y\ 


1  +  rj 


1+1 


d» 


14-80 


dx 


(16) 


ReDD 


1  +  n 

Integrating  this  equation  from  x  =  0  to  x,    and  noting  from 
equation  (11)  that  when  x  =  0 ,   >]  =  0 ,    we  obtain 
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nr 


•116 1\  +  16 


21| 


0  L    1 +  n 


(1+1  )/j/ 


%-2>7 

1+  1 


'14.-2  ^ 

16  Al + 


1+  y\ 


x+y\  1  +  >| 


5ol 


4-2  n 

1+  n 


1  +  A     _ 


d7| 


l|8o 


Diu 


D     '0 


dx      .  (17) 


The  integrations  on  the  left-hand  side  of  the  equation  above  are 
carried  out  term  by  term  with  the  aid  of  integration  tables 
as  follows. 


V  -116>1     +   16 


1  +  r\ 


24 


an 


'  132 


-ll6(r[  +1)    +  132 


/o        1 +  n 
A 


&n 


d>1     -  /       ll6dn  =   -116^     +   132   ln(l+>7) 
0  1+1  k  (18) 

d7] 


10  j  h-m 

(1+1) 


&y\   =  2k  {  — -? 


l+?j)  (I4.-2  y| ) 


1+  Y[ 

2k 


d  1 


0       /l/-2/|2+    2)|+i; 


=   2^ 


1  -lw+2 

—    sin"1    (— ^—) 
Y!  6      . 


-12  V  2 


.,-„-! 


-2»l+l - 

( ) 

3      Jo 


=  12  </? 


sin"1    ( ) 

3       Jr 


2i|d77  '  2)1-1  1 

\  =   12Y2  sin"1    (— — )    +   I2V2   sin"1    (-)     . 

0  A-21  3  3 

(i+7|)y (19) 

1+  >] 
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A  A-2JJ  ,        >■ 

-164/ d>i  =  -16    4/(4-2  ?i)(i+  n) 


+  3 


-16 


-16 


-1  -i+n+2 

sin"1    ( ) 

\fl  6 


2/j-l 


'l/l|+2'2-27l2   +  —   sin"1    (— — ) 
'  V2  3 


tl 


'V4+2?l-2>/   +  -p.   sin"1    ( )    -    2 


3 

?2 


3 


+  -_   sin"1    (-) 


1 
3  J 


k-2*) 


0       ri+ 1 


d>i     =32-24  V2   sin-1    (-) 

3 


V4+ 


2/1-1 


16  y  11+2*1-2^    -    2.1;    V2    sin-1    ( ) 


4-2?? 


1+  »                       A     -]/ 4-2  71 
—  d>j  =   2  T7?   d 

1   +   f\ 


=  4 

=  lj. 
=  k 


b      <1+^)" 

'     4/4-21  1  -4,,+2 

.  +  ^  sin"1  ( ) 

Vl+  >}  V2                      6      J  0 


-I/4-2I 
A/i+>\        42 
4/4-21         1 


1  2'J-l"1' 

-=   sin-        * 

^  3      Jn 


a-i 

_  sin"1    ( ) 

V1+  »1  V2  3 


1  .      1 

+  2  -—p.  sin-1    (— ) 
V  2  3  J 


(20) 
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1+  n 


i  +  n 


dn     =8-2   V~2   sin"1 


%-2Tr\  271-1 

l|f -    2    1/2    sin"1    ( ) 

1+U  3 


50  nf 


Let 

then 
50 


l(.-2>7 
1+  >-; 


1+1 


dl    =  50 


'1   ti-A 


^yi4.-2>7 


o    d+1) 


372 


d/i 


1   +   >1     =   v,  dv   =   d>]  ,  >1   =   v-l 


YlYk-^n  r    (v-1)  -j/\-2v+2 

d/'l     =   50   )    — —        -  dv 


(l+)?)3/2 


v3/2 


50  ( 


V6-2v  »   ^/6-2v 


dv   - 


=  50 


1/7  /     v3/2 

f -|/v(6-2v)    +   3  -/=   sin"1 


dv) 


(21) 


-kv+6 

( ) 

6 


I  -  V6-2v  1  -kv+6 

■2  |. ==—  +  -7=  sin"1    ( ) 

L        V7  fi  6 


,  / 3  i      1"2^ 

50     1/(1+1)  (L-^)    -  -=  sin"1    ( 

L  72  3 


+   2 


k-2*l 


i+  »\    yj 


L  i— -) 


1-21] 

3 


Thus 


5o  i    tjx  an 


50l/l4.+21-2n2  +   100 


1+  n 


2>1-1 

+   125^2    sin"1    ( ) 

■    •  3 


-fl 


5oVl|+2''(-2'l2   +   ioo  y  • 


1+  '1 


+   125 V2   sin"1    ( ) 

3 

1 
-    100-200   +   125V2    sin-1    (-) 

3 
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1^-2*1 


1    5°n1 
/-<  '   1+  '1 

Jo        1  +  >? 


d>1     =   50Yij.+2^-2n2   +   100 


1|-2  1 

i+  n 


+   12$V2    sin"1    ( ) 

3 

-   300   +   12572   sin"1    (- 


3 


(22) 


Now  substituting  equations  (18)  through  (22)  into  equation  (17), 
we  obtain 

-ll6n  +  132  In  (1+  Y[)    +  (12  Vi"  -  2I4.-/2  -  2  f2  +  I25V2) 

217-1 

sin"1  (— — ) 
3 

+    (12/2    -    24V2   -    2 1/2   +   12S-V£)    sin"1    (-)    +    (-16+50)1/ 14.+2>1-2>12 

3 


'14.-2*1 

+  (-i^+iooy +  32+8-300 

1+  1 


I460 


DRer 
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1+80  2)1-1 

X=   -1161     +   132    In    (1+n)    +   1111/2    sin"1    ( ) 

DReD  3 


1111/2    sin"1    (-)    +   3k-Vk+2>\-2r\2  + 
3 


96 -y 26 

1+1 


*  1  •  r  ,    2>i-i 

= -116  n    +   132    In    (1+n)    +   111V2    sin"1    ( ) 

DReD  I4.8O  L  3 


+   111  I/2    sin"1    (-)    +   341/V21-2H2  +  961/  -    260 

3  1+  1 


1 


16  /J4.-2I 


58      22  17   . 

Y\ In  {1+n) Vij.+2n-2l2 

DReD   16  [l5     5  15  5  '  1+  *? 


37  Y  2        2*1-1    37  V2        1    26 

sin"1  ( )  -  sin"1  (-)  +  — 

10  3       10         3  3 


(23) 


Thus  the  velocity  distribution  in  the  entrance  region  is  com- 
pletely defined  by  equations  (1),  (12), and  (23). 

The  Length  of  the  Entrance  Region 


When  the  flow  is  fully  developed,-;-  =2,  or  f\   =  1.   There- 

u 
fore  from  equation  (23)  we  obtain 

1 
f(1)  =  —  (3.86667  -  3.0ii.964.  -  2.26667  -  3.2  -  1.77776 
16 

0.i4.6li).l 

-  1.77776  +  8.66667)  =  =  0.0288  . 


16 


Thus  the  entrance  length  x   is 


0.02{ 


DRer 
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xe  =  0.0288  DReD 


(2k) 


Pressure  Drop 


Prom  Bernoulli  equation,  the  difference  in  pressure  between 
any  two  locations  in  the  entrance  region  is 


P]_  -  p2  =  -  Cu0)2  -  -  (u0)1 


(25) 


where  p^  and  p2  are  the  pressure  at  X]_  and  x2,  respectively, 
and  (uq)-l  and  (uq)2  are  the  velocities  at  the  central  core  at 
x-j_  and  x2,  respectively. 

From  equation  (11),  we  have 

(u0)1  =  ud+  n±)  ,        (u0)2  =  u(i+  n2)    . 

Hence  the  equation  (25)  becomes 


P-2 
Pi   -    P2   =  -  u 
2 

2 


d+>i2)2  -  (i+  n^2 


(>122  -  r^2)  +  2(n2  -  V 


(26) 


The  pressure  drop  at  the  entrance  region  and  at  the  fully 
developed  region  can  be  obtained  as  follows. 

For  x<  0.0288  DReD,  the  pressure  difference  between  the 
inlet  edge  and  any  location  in  the  entrance  region  is  obtained 
from  equation  (26),  as 
fu"2 


P0  "  P  = 


(>r  -  n02)  +  acvi  -  n0) 


Since  Y\. 
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a-2 

" u     o 

P0  -  p  =  (1  2  +  2»)   . 

u        2 

If  x  =  xe  =  0.0288  DReD,  7]  =  1.   Therefore  the  total  pressure 

drop  in  the  entrance  region  is 

3 


Ap=-  P 
2 


-2 
u 


For  x  >  0.0288  DRen,  the  pressure  gradient  in  the  fully 
developed  region  is  constant  and  is  given  by 
dp   32(j.u   6J+  fu2      1 


( ) 


dx    D2     ReD    2    D 
The  pressure  difference  between  the  inlet  edge  and  a  loca- 
tion in  the  fully  developed  region  is  the  sum  of  the  pressure 
drop  in  the  entrance  region  and  the  pressure  drop  between 
x  =  x„  and  the  location  considered.   Therefore 

3   2   6^   f>u2  1  /x 
p0  -  p  =  -  fvT   +  ( )  -     dx 

2  ReD    2    D  ^ 

3  -    bk        fu2      1 

=  -  Pu     +  ( )  -  (x  -  0.0288  ReDD) 

2       ReD    2    D 


P0  -  P        (>k 

U  _9   =  3  +  x  -  6[j.(0.0288) 

pu  ReDD 


p0  -  p     6J+ 

——2-  =  x  +  1.16   .  (27) 

P  u     ReDD 

2 

If  the  fluid  enters  the  pipe  from  a  region  where  the  pressure 
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ls  P<?o  snd  the  velocity  is  negligible,  there  is  a  pressure  drop 

f  u2 

of   between  that  region  and  the  inlet  edge.   Then  equation 

2 

(27)  becomes 

P_  -  P   *>k   x 


P u2      DReD 


+  2.16   .  (28) 
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APPENDIX  II.   LANGHAAR'S  LINEARIZING 
APPROXIMATION  METHOD 


The  Family  of  Velocity  Profiles  in 
a  Circular  Tube 


The  axial  direction  (x-direction,  see  Fig.  1)  component  of 
the  Navier-Stokes  equation  in  rectangular  coordinates  is 

0\x  0U  ?u      3\x  1  <3P      <32u    22u    3   u 

—  +  u  —  +  v  —  +  w  —  =  g--  —  +  y  ( — -  +  — -  +  — p 

St  3*  37  3z  ^       I    <?x     3^        3ld        3* 

(1) 

where  g   is  the  x  component  of  the  acceleration  of  gravity. 

(Reference  1. ) 

It  is  assumed  that  the  flow  is  steady,  incompressible  and 

laminar,  the  fluid  properties  are  constant;  the  transverse 

velocities  v  and  w  are  small  compared  to  u;  hence  the  pressure 

3V 
gradient  -  —  is  a  function  of  x  only.   Moreover,  the  velocity 

3* 

at  the  entrance  is  assumed  to  be  uniform  and  equal  to  the  average 

velocity  u,  and  the  term  — -  may  be  dropped  in  equation  (1)  in 

22u  <?2u 

comparison  with  and  — - 

<?z2     37 

With  the  above  assumptions,  equation  (1)  reduces  to 

3u  <?u     pu  1  dp     32u        a2u 

u  —  +  v  —  +  w  —  =  g   -  -  —  +  »(— -  +  — )  .    (2) 
3x     37  dz  r    dx     37l        32 

In  order  to  linearize  the  differential  equation  the  convective 

terms  on  the  left-hand  side  of  equation  (2)  are  replaced  by  the 

following  term: 
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,?u  #u  ?u  „ 

u  • —  +v  —  +  w  —  =  \>p  u  (3) 

<i>x  ^y  <?z 

where  p  is  a  function  of  x  only.   Equation  (2)  then  becomes  a 

linear  differential  equation, 

1  dp     #2u   J2u 
1)  (32u  =  gx  -  -  —  +  S>  (—  +  — -)  • 
f   dx     ^y2   ^z2 

Dividing  both  sides  by  V  yields 

1       1  dp     32u   a2u 
p2u  =  -  (g )  +  ( +  — -) 

v     x     />  dx        di2     az2 


a  2u   a  2u 

+ Q^u  =  -<  (1|) 

9       9 

?y^  ,?z^ 


where 

1       1  dp 

a  function  of  x  alone. 

The  assumption  expressed  by  equation  (3)  can  be  justified 
partially  on  the  following  theoretical  grounds: 

1.  Equation  (3)  is  satisfied  by  (3  =  0  in  the  fully  de- 
veloped region. 

2.  Equation  (3)  is  satisfied  at  all  points  on  the  wall 
irrespective  of  the  value  of  p  because  u  and  the  left- 
hand  side  of  equation  (3)  both  vanish  there. 

3.  If  it  is  assumed  that  the  laminar  state  is  attained  by 
the  development  of  a  boundary  layer  in  the  entrance 
region,  then  equation  (3)  is  satisfied  at  all  points 
in  the  central  core  of  the  fluid.   This  conclusion  is 
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reached  from  the  fact  that  —  and  —  in  equation  (3) 
vanish  over  the  flat  central  portions  of  the  velocity 
profiles.   Accordingly,  equation  (3)  becomes 

u  —  =  Vp^u 
dx 


or 

3u 


-=VP^  .  (5) 

3x 

■  au 

Since  u  and  —  are  functions  of  x  alone  for  the  poten- 
tial  flow  in  the  central  core  region,  it  follows  that 
equation  (3)  is  satisfied  by  equation  (5)  in  that  region. 
k-   Equation  (3)  is  satisfied  at  all  points  of  the  inlet 
section  x  =  0,  since  the  thickness  of  the  boundary 
layer  is  zero. 
5.  Equation  (3)  is  justified  except  for  those  points  which 
lie  within  the  boundary  layer.   In  the  latter  part  of 
the  entrance  region,  a  distinct  boundary  layer  does  not 
exist  since  the  central  parts  of  the  velocity  profiles 
are  appreciably  rounded.   However,  in  that  region  the 
flow  approaches  the  fully  developed  laminar  state  for 
which  equation  (3)  is  valid,  as  explained  in  1. 

In  cylindrical  coordinates,  we  have 

poo 
y^  +  z-  =  r* 

u  =  u(r) 

au   du  }v        du  y 

jy        dr  Jy        dr  r 
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32u   d       du  y    y  3       &v.    av       du  9      y 

=  _  ( )  = (— )  —  + (-) 

31  ay   dr  r    r  3T      dr  ay   dr  <?y  r 


r  -  y  — 

y2  d2u    du  3  J 

9     ?  ? 

r  dr    dr  r 

y2 

?   2        r  -  — 

y11  d  u    du  r 

~2  Z5  + 


r^  dr    dr 

Similarly, 


r  -  — 
53  u   z   d  u   du     r 


^z2    r2  dr2    dr   r2 


Therefore 

32u       32u.       d2u  y2  +  z2    du  1       y2       z 

gy2   gz2        dr2    r2       dr  r2      r        r 

d2u    du  1         y2  +  z2 

=  _  +  -  _  (2r ) 

dr^    dr  r^  r 

d2u   1  du 

dr2    r  dr 

Equation  (Ij.)  then  becomes 

d2u   1  du 

+ B^u  =  -<  .  (6) 

o 
dr11   r  dr 

The  particular  solution,  by  inspection,  is 

u  =  -  -5  . 

P2 

The  complementary  equation  is 


(7) 
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d"u   1  du 

+ p2u 


(8) 


dr    r  dr 

which  is  a  modified  Bessel'a  equation  of  zero  order.   The  solu- 
tion of  equation  (8)  is  (Ref.  3) 
U  =  Cx  I0(pr)  +  C2  K0(pr) 
where  IQ  is  a  modified  Bessel  function  of  first  kind  and  of  zero 
order,  and  K0  is  a  modified  Bessel  function  of  second  kind  and 
of  zero  order. 

Since  KQ((3r)  approaches  to  infinity  while  u  remains  finite 
as  r — >  0,  therefore  C2  must  be  zero.   Thus 

u  =  Cx   I0(pr).  (9) 

The  solution  of  equation  (6)  is  the  sum  of  the  particular  solu- 
tion (7)  and  the  complementary  solution  (9).   Hence 


ci  W')  -  — 


(10) 


From  the  boundary  condition  that  u  =  0  at  r  =  r0,  we  have 
0  =  01    I0(Pr0)  -  1 


or 


where 


~2    '   °1  X0(y) 

y  =  Pr0  . 
The  continuity  equation  is  used  to  eliminate 


(11) 


1        1  dp 
-  (SX  "  -r  — ) 

n     x    C  dX  J 


The  continuity  equation  is 


kpq      u  =     j  2nrudr 

'0 
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2   - 

-  rn      u 


^0 


rudr 


Substituting   equation    (10)    into   this   equation  yields 


-  rn      u 


<r0 


cl   Io(Pr>    -  —    jd 


/-r0  -<        fr0 

C^     I         rl0(pr)dr   -  -x  rdr 

'0  P       'o 


1  0      f"  .      .'o 


(pr)    I0(pr)d(pr) 


rdr 


(pr)    In(pr) 


L 


if)    1^?) 


0         P2         2 

,      „   2 


P2      2 


P2r02   u  =  Cx 


C/)    I^V) 


pV 


or 

o 

-y2  u  =  c  y  i  (y)  -  _  .  i-    . 

2  X  p2         2 

Substituting   equation    (11)    into   equation    (12),    we    obtain 


(12) 


-f2  u  =  CjVl^/J    -Lca   I0(/) 


8o 


Dividing  both  aides  of  this  equation  by  //2  yields 

°l[?llM  -  ioH  -S  . 

Using   the   relationship    (Ref„    Ij.) 

i2(y)  =  io(V)   -  -.Ii(V) 

7 


/ 


I-L^)    -    I0(/)    =    -    I2(/) 


we    obtain 

-   C1I2(/)    =  u 


Cl   = 


Substituting  equation  (13)  into  equation  (11)  yields 


InO') 


(13) 


(Hj.) 


^         i2m 

Substituting  equation  (13)  and  (ll|.)  into  equation  (10),  we  ob- 
tain the  equation  for  the  velocity  profile  given  by 

loM  -  io(Pr)  _ 

U  =  U  -  u  . 

1 2  00      X2(^) 

Dividing  both  sides  by  the  average  velocity  u  yields 
u  _  I0(/)  -  Iptpr) 
u      i2(Y) 
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i0(r)  -  i0<r— ) 

u  r0 

_  = .  (15) 

u  i2(r) 

Thus  if  the  relation  between  /  (=  prQ)  and  x  is  known,  equation 
(15)  completely  defines   the  velocity  profiles  in  the  entrance 
region.   The  next  task  is  to  find  the  function  Y (x) . 

The  Function  /(x) 

The  momentum  equation  can  be  obtained  by  integrating  the 
Navier-Stokes  equation  over  the  cross  section. 

In  the  cylindrical  coordinate  system,  the  simplified 

Navier-Stokes  equation,  equation  (2)  becomes  (Ref.  2) 

2u     gxk  ■       1  dp     tf2\x        1  ju 

—  +  u  —  =  gx-_  —  +V>(— -  +  -  -  ' 
<?r     ^x        f  dx     Jr"5   7  01 

and  the  continuity  equation  becomes 


vr  —  +  u—  =  gx----+V>  ( — -  + )  (16) 


3vr   vr   gu 
3r  r    #x 


1  0         0u 

(rv  )  +  —  =  o  (17) 

r  ^r        ^x 

where  vr  is  the  velocity  in  the  radial  direction. 

Integrating  equation  (16)  over  the  cross  section,  we.  have 

fr0  <5u     ju 

'  rr   —  +  u 
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r°  /r0   1   dp  ,r0     d2u       1  0n 

gx27irdr   -   / 27trdr  -t-  /      V( + ) 2nrdr 

0        C   dx  /0  jr2        r  2v 

(v0  $U  ?u  (v0  i   dp      /-r0 

2  J        (vr  —  +   u  — )rdr  =   gx  I       d(r2)--—  d(r2) 

0  <?r  ^x  JQ  fell    /0 

,    /r0      _j2   u        x    ?u 

+   2U  I  ( + )rdr 

'0  3  ?2       v  fr 


rTQ 


(vr  —  +  u  — )rdr  =   g  r. 


2        1   dp        2 


'0  3v  3*.  f  dx 

/^O    '<J2u        !  gyx 

+   2»  I  ( + )rdr      .         (18) 

)q  d?2        r  Jr 

The   integration   of   the   left-hand   side    of   equation    (18)    can   be 

simplified   as   follows: 

r-    .t* 

3u 


'r°.     ^u       ^u  r  r1"0    ^u  rr°    ^ 


2  I        (v     —  +  u  — )rdr  =   2 
!0  2?  3* 


u  —  rdr 

/-ro       5u  1   fr0   ^(u2) 

I        vr  —  rdr  +  -  /         rdr 

0  9T 


2/0  dyi 


Let 


then 


(19) 


5u 
7? 


du'    =   d(rv   )    and   v'    =  u    . 
Thus   integrating  by  parts   yields 


'0        3u 

vr  —  rdr   = 
3? 


■0         ,p0 
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[0  3u 

Jr 
From   equation    (17),    we    obtain 

r  3r  fix 

Hence   equation    (20)    becomes 


rLQ        1   3(rvr) 

u rdr 

0  r       <?r 


.  '0  2u 

I  vr  —  rdr 

'0  ^r 


*u 


)       u  — -  rdr  =  - 
0  3*  2 


'r°   ^(u2) 


rdr 


'o         5x 
Substituting   equation    (21)    into   equation    (19)    yields 

2   I        (v 
'0 


(20) 


(21) 


3u  0m. 

•r  —  +   u  — )rdr  =   2 


P°^(u2)  i    rr0  5(u2) 


.2    '0         dx 


rdr  +  — 
2 


ax 


rdr 


=   2 


rdr 


;0  <?x 


3U  2u  d 

(vr  —  +  u  — )rdr   =2  — 
0  r  ,?  x  dx 


p 

u   rdr 


Thus    equation    (18)    becomes 


dx 


2  —  I       u^rdr   =   g  r( 


xx0 


1   dp 

}  Tx 


r0  32u 


(22) 


1   2U 

+     2)>  (__    + ); 


0        <?r 


r  <5r 


Prom  equation    (15),    we   obtain 
3u 
3r 


U         0  f 

— I0(—  r) 

I2(r)   <?r     U    rQ 


u  9/  / 

(— )    I-,(—  r) 

I2^)      r0  r0 


ai*. 


2u 


up 


dv  I2(Y)      1 


In(Pr) 


and 


a  <^u 


3** 


up2 
I2(/) 


I1*(PP) 


(23) 


(2^) 


The 


refore    the   last    term  of  equation    (18)    becomes 


fv0    S2u        1   2u  v  fV0 

2Vj        ( + )rdr  =   -21) 


Iq       dr         r  ?r 


up2 

ioC/) 


Ij    (pr)    rdr 


-2V  I,  (pr)dr 

J0    i2(r) 


y  _     v 

U  (■  |  u  f 

=  -2V*  /     (pr)I-,    (pr)d(pr)  -2^ I-,  (pr)d(pr) 

i2(/)/0  l2(f)l0 


L2vr  I  /Q 
2Vu 


IoC/) 


(pr)In(pr) 


,1  | 

I1(pr)d(pr)    +       I1(pr)d(pr>f 
0  ' 


,,     /-r0     <?2u          1  ju 
2V     I  (— -)    + )rdr 


J0       'o 
2^u 


/i^/) 


(25) 
/0  ^  r  <?r  I2(7) 

Substituting   equations    (22)    and    (25)    into   equation    (l8), 

we   obtain   the  momentum  integral   equation 


2liu/l1(-/) 


+   Tn*    ( 


0       ^x 


1    dp  d 

— )    =    2   - 


c 


p 

u   rdr 


dx 


dx 


Since  the  radial  velocity  is  zero  in  the  central  core,  eque 
tion  (16)  becomes 


(26) 


1  dp 


du 


,?2U    1  JU 


Z  —  =   _£  -j  f      +1—1) 


dx 


r"        r  3?   r=0 


(27) 
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Prom  equations  (23)  and  ( 21(.)  ,  we  have 
1  3u      up    Ii((3r) 


r  dv  l2{i) 


T 


lira      t1   3u,       UP   iim    '^CPf)  UP2  iim 

r  ^r      l2(/)         P  I2(y)r^U   ! 

u(32     ,         up2  I0(0)  -  I2(0) 
II  (0) 


I2C)   X      _    I2^) 


1   up2 


and 


2   I2(0 


UP2     ,,         up2   I„(0)  -  IP(0) 


(28) 


Ii  (0) 


r2   r_°     I2m      l  I2(/) 

1   up2 


(29) 


2   I2(V) 

Because  of  equations  (28)  and  (29),  equation  (27)  can  be 

written  as 

1  dp      du0    Vup2 

g =  un  +  

x   P   dx    u  dx    !„(•/) 


o 

1  dp      duQ     V  u/. 

f  dx      dx    rQ2I2(/) 
Substituting  this  equation  into  equation  (26),  we  obtain 
2^/1^/)         duQ    yuY2     d    fr0 

777 +  r0  u0  +  =  2—1    u2rdr 

^^  dx    I2(y)      dx  /- 


»u/2   -    2»u/I1(/)  d        rP° 

0 
/r0 


Ipoo 


1  d 
u2rdr   -   rQ2  (uQ2) 

2  dx 


=2  —     I  u^rdr   

dx   ''g  dx 


2 
Ug   rdr 


dx 


(2u2    -   uQ2)rdr 


i2  -  2yi1(y')        i    d 


1 2  00 


V  u   dx 
urQ2    d         f1 


(2u2   -   u02)rdr 


<2  Z?  -  Z1-X- )d(  — )       .  (3D 

V        dx    IQ  u  u  r0        r 


Let 


u  u0  r  urv  x 

U   =  -   ,   Ug   =  — ,    H  =  —   , ,  Rer  =  — -  and  X  =  

u  u  rn  ))  i"nRe, 


Then  equation    (3D    becomes 


Y2    -   2-/l1(-/) 
Io(/) 


=   r0   Rer 


dX 


(2TJ2    -  Ug2)RdR 


J   dx 


—     /       (2U2   -  Un2)RdR 
dX    /0 


2YI1(<y)    -  y2        d        f1 


!?(/) 


-  =  —  (Unc    -    2IT)RdR 

dX 


Noting   that   equation    (l£)    can   be   written   aa 

Iq(/)   -  io(/R) 


U   = 


lo(f) 


(32) 


(33) 
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the  integration  on  the  right-hand  side  of  equation  (32)  can 
be  performed  as  follows: 


(Un      -    2lT)RdR  = 


i0(r)   -  I0(0) 


Ic2(/)  2   L      Jn 


Ip(/) 


loC)  -  io(/R) 


Ip(/) 


RdR 


[l<>V) 


o      i22(/)  £ 

2I0(/)I0(/R)    +   I02(fR) 
ri02(/) 


RdR 


I22(/) 

[ip(/)   -  i]2 

2I22(f)  I22(/)L      2 

-    2I0(/)    —  j        (/R)I0(/R)d(/R)    +     |        IQ2(/R)RdR 

7    'o  '0 

[ioci  -  ij2    i02^)    i^ott')  r        f 

2x22(/)  -  -^ + /-%^  h>ii^>jo 

a 

I02(fR)RdR 


l22^)     ^0 

EW     -    l]2        iQ2^)         iHO(?0    Il(f)  2 

2I22(/)         "    I22(r)    +       /I22(/)  '"    I22(n 

[       I02(/R)RdR 


(3W 


For  the  integration  of  the  last  term  of  equation  (3k),    the 
following  formula  is  used  (Ref.  3). 

.2 


/   xl„  (Xx)dx  =  - 

'0  2  l 


I'  (Xx)  -  (1  +  — — )  I  2Ux)   .   (35) 


2    2    n 


Thus 
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Iq    (/R)RdR  =   - 


R"  .     .  o  o  -i  1 


i  2 


(i0'Vr)  -  (i  +  — _)  i02(yR)} 


y^R 


Jo 


-;  [V^  -  io2w] 


I    I02(/R)RdR  =  -    [l02(/)    -   I0'V)]  =  -    ["l02(/)    "   Ii2(/)]    • 
/n  2  2    ^~ 


Equation   (3i|)    now  becomes 


(U02    -   2U2)RdR 


[l0-l]2        I02        Ijloll  2  1 


21  c 


/I; 


(36) 


do2-!!2) 


tln-D' 


2Ir 


?I22  I22        I22       °      2 


^o2-ii2 


I22 


(lQ-D: 
2I22 


2InI 


0X2 


I?2 


(U0^    -    2U^)RdR   = 


(I0-1)2    -   lj.I0I2    +    21] 
2 


'0  2I2- 

in  which   the'  argument   of  I 's    is    y  . 

Substituting   equation    (37)    into   equation    (32)    gives 
2/1 i   -  /2        a        (I0-1)2   -   ^I0I2  +  2IX2 


(37) 


Let 


dX 


2I22 


f(/) 
g(/) 


1|I0I2    -    (Iq-1)2    -    2I12 
2I52 


2/ln    -  /2 


(38) 

(39a) 
(39b) 
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then   equation    (38)    becomes 
1  d 

g(-/)  dX 


1               d     r        n 
_  = |f(/)j 


g(/)  d/ 


[f(y)] 


d/ 


dX 


i  dy 
=  -  fi(y)  _ 

;(/)  dX 


(39) 


This  is  a  differential  equation  which  expresses  f    as  a  func- 
tion of  X. 

For  a  large  value  of  Y  ,    1    [l/)    can  be  represented  by  an 
asymptotic  series,  as  (Ref.  l\.) 


In»0 


lj.n2-l2    (ij.n2-l)(^n2-32) 


r  ! + 

(2n/)1/2  L    8/-1I        (8/)2  2! 

I7i[1+^(^J/ 


(2rcr) 

where       £•„(/)    =  0    as    / — >  oo  . 
Accordingly,    from   equation   133) ,    we   have 

I0(V)    -   I0(/R) 


lira        TI    _        iim 
/ >=^  / »  o& 


I2('/) 


Am  Am        X0(/R) 

/R 


/-^°°I2(/)        ^°*    I2(/) 


a.               (2TI/R)1/2 
t  *  Im       


(1*0) 
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,         /im 


V5  y^°"  T^rT     ' 


e 
Since 

0  <  B  5  1  or  1  -  R  >  0 

^?-  u  ■ x  -  jf  °  -  x  ■ 

But 

u 
U  =  -  =  1 

u 

at  the  inlet  section,  x  =  0.   Therefore, 

r(H^    X(/)  =  0  . 
Thus  integrating  equation  (39)  from  X  =  0  to  X  leads  to 

-x       y 

dX  =  -  I   g(/)  f '  (V)6Y 

'0  L> 

or 

X  =  f     g(/)  df(/)   .  (lj.1) 

V 

For  the  numerical  integration  of  equation  (Ip.),  the  following 
change  of  variables  is  made. 

Z  =  £(/)   and  Y  =  g(/)  (lj.2) 

The  values  of  Z  and  Y  as  / — >  °o   can  be  found  as  follows: 
From  equation  (1+0),  we  know  that  as  /->»o 
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These  relations  will  be  used  in  the  evaluation  of  the  values  of 
Z  and  Y  as  '/ — i  o^>   as  follows: 

o.  -       ,  f.  ^0I2  "  do"1)2  "  2Il2 

1  iim   2I0-X   X 

=  2 i  +    .fim     — —  =  -  , 

2  y-^*>    2I22     2  ' 
because 

1 

Am       f^  _        Am       !_j£  _    ^im       i^  _  . 
r-^c^2l22         *-»»*      l22    -?-*-»     2l2 


?im       „  _       ./im        _u/\    .        Am 


Y=        *     .        g(/)    = 


r — ><>=>  / — 5b»  &  y—^o^  2/1-,-/2 

Am       _       Am 


■/—***      Xl     ^      ^-*o^ 


^2 
2  7—   -   —  2/ 

I2        12  12 


Since 


A»       2l-        Am       _il__    f2_,l/2        /im       zf_  a  0 


(2ti/)1/2 
as   L'Hospital's   rule    is   applied  repeatedly,    therefore 

>im  y  =       -Am        _  q 

Next,  we  shall  find  the  limit  of  Z  and  Y  as  /— >■  0.   I  functions 

may  be  expressed  in  power  series  as  (Ref.  i|) 


92 


In('/)  =  i"n  Jn(i/) 


(  i   /)n  f    (  -L  /)2  (  ±  ^)4 


1  + 


i!       I 


l(n+l)     1  •  2(n+l)(n+2) 


+  • 


(^-/)6 


1  •  2  •  3(n+D(n+2)(n+3) 


Let 


then 


y  = 


In(/)  =  1  +  a2  +  —  +  —  + 
1+    36 


I-iOO 


I2(/)  =  — 
2 


a2   a*J-   a6 

1  +  —  +  —  +  +  . 

2     12    XJ4J4. 

a2   ak        a6 

1  +  —  +  —  +  + 

3    24   360 


Thus 


lj.I0I2  =  2a2  (1  +  -  a2  +  -  ak  + 

3     8 


(I0  -  l)2  =  ek   +  —  + 
2 


21-,^  =  2a' 


1  +  a2  +  —  a^  + 
12 


2   a^  r    2  9 

21  2  =  _  h  +  _  a2  +  _  _  _ 
2   L     3 
Now  Z  can  be  written  as 

z  _  4I0I2  -  (Iq-1)2  -  2IX2 
2I92 
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P   8a^   5b6  .   a6  5   , 

(2a^+  +  +...)-(a^+  —  +...)-(2a2+2a^+  -  a6...) 

3  k  2  6 

7k       ^ 

2^3 
a^   a6 
3    12 


a4    a° 

—  +  —  +  .  .  . 

2     3 

Since  /  >  0,  as  a >0,  we  have 

a^   a6 


y — >o 


z  = 


a »0 


3    12 


a^    a6 

—  +  —  +  ... 

2     3 

1    a2 


3   12 


Similarly  we  have  for  Y 


Y  = 


12 


2        3 


—     1  +  — •  +  —  +    . 
2      i         3  2k 


2Ylx-i2 


a2 

2  6 


83  5         5 

+  —  +  —  a->   + 

2  12 
9^        H6 


-   k&< 


2a^  +  _   a6    +    .., 


and 


9^ 


A   «6 


,iim   v  _  Am 
/— >  0  Y  ~  a >  0 


2    6    lj.8 

( )  =°o       . 

+  — 
6 


2aV-  '  -6 


Summarizing,  we  have  the  following  limits  for  new  variables 
Z  and  Y: 

Z >-   ,         Y- — >  0   as   / >  cxfl 

2 

2 

Z  — > ,    X-^>^o   as  / >  0   . 

3 

Thus,  equation  (i(J.)  may  be  expressed  as 

rl/2 

X  =  j     YdZ  .  (ljj) 

Equation  (lj.2)  may  be  regarded  as  a  parametric  representation  of 

a  curve,  and  equation  (I4.3)  shows  that  X(V)  is  the  area  under 

1 

this  curve  in  Y-Z  coordinates  in  the  interval  (Z,  — ) .   Therefore 

2 

if  we  assign  a  set  of  values  to  /  ,  a  table  of  points  on  this 

curve  may  be  calculated.   Then  by  numerically  integrating 
equation  (1|3),  we  obtain  the  corresponding  value  of  X. 

In  this  manner,  values  of  X  corresponding  to  various  values 
of  "/   are  calculated  and  listed  in  Table  1  (Ref.  1).   Thus,  by 
giving  a  value  of  /  ,  X  can  be  found  from  this  table.   X  being 
given,  the  velocity  profile  of  this  corresponding  section  is 
then  defined  by  equation  (15)  or  (33)  •   Hence  the  velocity  field 
in  the  entrance  region  is  completely  defined. 
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The  Length  of  Entrance  Region 

As  usual,  the  length  of  the  entrance  region  is  defined  as 
the  distance  from  the  entry  to  a  section  where  the  center  line 
velocity  is  99  per  cent  of  the  fully  developed  one;  that  is, 
u0  =  0.99  (u0)F.Di 

uo      (uo)f.d. 
—  =  0.99  . 


(uq)f.d.  __ 


However, 


hence 

u0 
U0  =  -X  =  0.99  x  2  =  1.98. 

u 
This  is  equivalent  to  saying  that  the  entrance  length  is  the 
distance  from  the  entry  to  a  section  where  the  dimensionless 
center  line  velocity,  UQ ,  is  equal  to  1.98.   The  values  of  UQ 
at  various  sections  can  be  obtained  from  equation  {33)    with 
R  =  0;  that  is, 

i0(V)  -   i 

0      I2</) 
Values  of  U0  are  listed  in  Table  1.   Prom  this  table 
Xe  =  0.227  when  UQ  =  1.98, 


Xe  =  =  0.227 

Rq  Rer 


=  0.227  , 


ur0 
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uD 


0.227  . 


Thus  we  can.  conclude  that 


D  ReD 


=  0.0^68   . 


The  Pressure  Drop 


dp 
In  the  entrance  region,  the  pressure  gradient  —  is  higher 

dx 
than  in  the  fully  developed  flow,  not  only  because  of  the  in- 
creased friction  loss  but  also  because  of  the  increase  in  kinetic 
energy  of  the  fluid  as  it  passes  downstream. 

The  method  which  is  employed  to  calculate  the  pressure  drop 
involves  an  energy  balance  in  the  entrance  region.   The  power 
which  is  dissipated  by  the  fluid  friction  in  the  interval  x  =  0 
to  x  =  x  is  given  by  the  dissipation  function,  as  follows 
(Ref.  2): 


E   =   |j. 


i' 


}w 


2  (— )   +  2(— r  +  2(— r  +  (- 


<?V 


?U 


3x 


+    (—  +  — 

di     3z 


)2+  ( 


31 

2U         g\i 


-y 


J-L 


'X 


dV 


3x        'pj 
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a.  I     dx 


.3  u  „     <?  v  „     ^w„ 

2(— r  +  2(— j2  +  2(— r 


ax 


^z 


(—  +  —)' 

^x   ^y 


+  (- 


d  +  (—  +  —  )2 


dA 


3H    i>U 
(—  +  — ) 
3  7       2z  3  z       J  x 

where  A  denotes  the  cross-sectional  area.   It  is  assumed  that 

w  and  v  may  be  neglected  in  this  equation.   Then,  by  the  con- 

,?u   ,?v   ^w  3'cl 

tinuity  equation,  —  +  —  +  —  =0,  we  conclude  that  —  =0. 
<?x  31  '  <?z  <?x 

Thus  the  preceding  equation  reduces  to 


E  =  \J.    j   dx 
'0 


,?U  o      ,Jli  , 
(  — )    +   (— ) 


dA  . 


(W 


During  the  time  interval  dt,  the  net  work  of  pressure  on 
the  body  of  fluid  in  the  interval  x  =  0  to  x  =  x  is  uA  (  Pq  -  p)dt. 
The  work  of  friction  on  this  body  is  -Edt,  where  E  is  given  by 
equation  ( J4J4.)  .   It  is  assumed  that  the  tube  is  horizontal,  so 
that  the  work  of  body  force  is  zero.   Since  the  work  of  all 
forces  is  equal  to  the  change  in  kinetic  energy,  we  have  the 
following  relation". 

'1 


uA(pQ  -  p)  dt  -  Edt  = 


f>\       (u3  -  u3)dA 


dt 


uA(p0  -  p) 


,x 

u.  1   dx 


A 


<?u 


<?u 


(— )    +  (— )' 

QY  3z 


1 
+  — 

2   'A 
In  cylindrical  coordinates,  we  have 


f  j     (u3  -  «3 ) 


dA 


dA 


(45) 
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gu      ?u  ,?u      <?u 

—  =  —  cos   & ,   —  =  —   s  in  9 
3  y       ,?r  ^  z       ^r 


Hence 

(— ) 
,3  y  ?  z  J  r 

Prom   equations    (23)    and    (35),    we    thus   obtain 


3  u   p 
(— )2   +    (— )2   =    (  — )       • 


2u 

-) 

'y 


(— )2    +    (— ) 


(— ) 
<?z 


2    R2         rr0 


dA   = 


u^    p 


*u 


5>u   2 


(— r  +  (— ) 

Al.3J  3z     . 

u2    p2 
=   2n 


dA    =   2k 


u2    p2    rr0 


I1  2  ( (3r)  2itrdr 


I12(pr)rdr 


L2        '0 


Io2(-/)L2 


■(-   I1'2(pr)  +  (1+  _ Jl^ffir)]     ' 


—2    2       2    - 

3tU    3> 


i22f)    L 


-V  (p*0)  + 


(3^r£ 


l  +  P2/p2 
P2-02 


ii  (P'o) 


l22(^) 

HU2 
Ip2('/) 


I22(^) 

Ttu2/2 

i22(y)  L 

Ttu2/2      r 
I22(>/)    ,. 


-{/ix(/)]2  +  (i+y2)  i12(/) 
-{-ix(/)  +  ^VO}2  +  (i+-^2)i12(/) 

/2I12(/)    +    2/I1(/)I0(/)    -   /2IQ2(/) 
2I1(/) 


1-fif)    +   IQ{?) 


/ 


I02(/) 


Ix2(/)     +    I0(7/)fl0('/)     -    12(/))      -    I02(f) 
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n2,/2 


7iur> 


Ip2(/) 


Xl  &)     -    IoM  I2(/) 


Thus, the  first  term  of  the  right-hand  side  of  equation  (I4.5) 
becomes  ' 


(J,    I       dx 
0 


2u    , 
I—)'    +    (- 


,?u 


i,A 


-2  rx  y2   r  2 

^U  7T    LJ1      "    I0I2 

-2    /^x    ^    r    ?  i         ^ro 

l*KU  —  [l^    -   I0I2]d(r0   —  X) 


=  /'nr   2  u^ 


.X      /2 


rs  f1!2  "  v2] 


Inlo      dX 


(46) 


'0        J-2 

The  second  term  of  the  right-hand  side  of  equation  (J4.5)  becomes 

1       (  -  P  -        fT° 

-f         (u^    -   U3)dA   =  -  u3  (tj3    -    l)27trdr 

2  ;a  2     ;0 

=  ti/^u3     I        (u3   -    1)  (r0R)d(r0R) 

,1 
=   7t,"u3r02     /       (u3    -   l)RdR      .  (^7) 

-  'O 

Substituting  equations  (I4.6)  and  (lj.7)  into  equation  (lj.5)  yields 


X  v2 


-o2(po  -  p)  =Ar02u3|    L-Jy  .  IoIj 


dX 


0    ij 


+  f%r02   u3     j       (u3    -    l)RdR      . 
0 


Thus 


/X 


/2 


P0-P   =  /V  j       — -   j^I12    .    I0I2|  dX  +  fu2  j        (u3-l)RdR 


100 


P0     -     P  „    fX    Y2         ^     n 


=2    — -  I I12  -  I0I2   dX  +  2  /   (u3  -  l)RdR 


I  ._„      /n   Io2  L  "     "  ^J 


0-2  '0   l2 

2 


7^ 


According  to  the  Bernoulli  equation,  the  pressure  drop  between 

a  reservoir  where  the  velocity  is  small  and  the  inlet  section 

is  —  P -a    .      Therefore , the  total  pressure  drop  between  the  supply 

2  1 

reservoir  and  section  x  is  -ft2  +  (pQ  -  p)  ,  or 


1  +  2  I   — -  I]_2  .  I0I2  dX  +  2  J   (U^  -  l)RdR  .   (24.8) 


;   t  2  ^ 

1/)-2  'o   I2   b  '0 

— pi 

2 
Let 

*(/)  =-  k2  -  i0i2 

i2 
then  we  have 

/-X   y2  ,X 

/    7T  L1^  -  V2  h*=     0M 

yo  .   i2  J     Iq 

The  evaluation  of  this  integral  is  difficult  due  to  the  fact 
that  <j>{V) — >oo  as  X — >0  (/-)oo),   Therefore  we  make  use  of 
a  limiting  relation.   Since 

klQl2   ~    do-D2  -  2Ii2 


dX  .  (49) 


ft/;  = 
ft/) 


2I22 
I22(l0l2)'-(loI2)(I22)' 

1  I22  •  2(lo-l)lo'-(Io-l)2(I22)'    I22(ll2)'.Il2(l22)' 

2  "  l2k  I? 
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I22(l0l2'+l2lo')-2lol22l2'         l22(Io-l)Io'-(Io-l)2I2l2' 


1^ 


X? 


-    2 


I22llll'    -    Il2I2l2' 


,         I]_  2 

Substituting  the  relations  I]_   =  Iq  -  — r  ,  1 2 '  =  Ii  -  -  I?  , 

7  / 

and  Iq   =  I]_,  we  have 

2  2 

i23f(/)  =  2  (i2[i0(ii  -  -  i2)  +  11I2]  -  2i0i2[ii  -  3  ia]J 

,  2 

-  I2(i0-Dli  +  do  -  2Io  +  Ddi  -  7  i2) 

Ji      9     2 

-  2I-|Ip(In  -  — )  +  2I12(I1  -  -  I,). 


Ll-L2^0 


1  Vil 
21, 


Simplifying  and  using  the  relation  =  Iq  -  I2,  we  obtain 

I23f(/)  =  -  6I0I1I2  +  3l!l22  +  I02Ii  +  2Ii3  +  y   I0I22 


"  y  J0   X2  +  Ixl2  -  2I01!  +  J?  ^^  +  Zl  "  7  X2  • 

(50) 
For  very  large  values  of/,  the  functions  of  I ' s  are  expressed 
in  asymptotic  series  as 
J 


i0(/) 


Y^j/ZkV 


I    2       /         9  2 

y    +  -  +  —  +  /  £ 

L  8       128 


o.w], 


=/ 


ix(/)-        . — 

y2 1/211/ 

e^ 
I?('/)- 


'  /2         3         l5  2 


r^-y^ 


12E 

15  105 


(5D 


8  128  2 


&m 


where     y_^  ^  /^n(/)    =0,        n  =  0,    1,    2. 
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Substituting  equation  (51)  into  equation  (£0)  yields 


9      15  105 

8  128 


f  (/) 


-6h*  +  -  +  — 


r-2    J      l5lr?    15      1051 

y    -  -  / y2 /  +  — 

L  8  I28jl-  8  128J 


+  3    /•  -  -  i 


'jz     v       9 

r    +  -  +  — 


15 

128 

2_ 


128. 


L  ft 


1* 

+  — 

y 

2 

"7 


'   2        / 

/      +  -  + 


128 


15  105 

8  128 

15  " 

12. 

15  105 


+    2 
2 


r2  -  -  Y — 

8  128. 


L  8  12" 


V 

8       128 


2      15  105 

8  128 


+   =T~  F(/)/ 


(52) 

where  F('/)  is  a  polynomial  function  of  Y  resulting  from  the 

substitution  of  equation  ($1)    into  the  expression 

k  2 

I]_l2  -  2IqI-,  +  —  Iq-^2  +  -'■I  "  v  ^"2  ^n  equation  (50)  . 


Since 


y — *<>o 


y      \J  ~  tL  ~  ±         y 

i2  l/iity 

y F(/)  =  0,  equation  (52)  yields 


9   ^     i°5  ■>  1,      , 

(/2 y  +  — r  f ■  (/)  =  y*  +  k'/3  +  .  .  .  , 

8      128 
after  simplification,  where  k'  is  a  numerical  constant.   Prom 
this  we  obtain 

f<(/)  =^  +  k^+  .  .  .  ;  (53) 

where  k  is  a  numerical  constant.   Moreover,  substituting 
equation  (51)  into 
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g(V) 


2ri,  -r 


2 


yields 

s(V)  =  —  +...  .  (5k) 

2/ 

Therefore,  from  equations  (J+l)  ,  (53),  and  (5Wj  we  obtain 

.00 

(/"3  +  .  .  .)d/=  — 2  +  .  .  .   . 
2  //  1^ 

It  follows  that 

/ — >=><=?  '  A  _  r 

*  2      i 

It  is  pointed  out  by  Langhaar  (1)  that  setting  f   X  =  —  in  the 

range  7  £.    100  is  a  good  approximation.   With  this  approximation 

the  integration  of  equation  (1|9)  becomes  easy  and  is  carried 

out  as  follows: 

.-X  JO. 000026         ,X 

|   0(/)dX  =  |         0(/)dX  +  J        0(/)dX 

'0  '0  '0.000026 

where  X  =  0.000026  corresponds  to  Y   =  100.   From  equation  (51)  , 
we  have 


1 2 

•1 

2 
+  —  + 

3 

• 

^1 
1? 

y^> 

1 

3 
2Y 

15 

+    o 

8/2 

+ 

lOlj. 


Hance 


3       6 


II2       l0       !       3         5 


i22      i2      y      y2     /3 

Making  the  approximation  that  y-X  =  —  for/  >  100,  we  obtain 
'    2  ^ 


— '/ +  3  +  - 
/ 

•~  -  X"1' 2  +  3  +  iox1'2 
2 


and 


^0.000026  ,0.000026 

i^(y)dX  =  /         (-  X"1/2  +  3  +  ioX+1/2)dX 


0  '0 

=  0.00518. 
Therefore, 

<x  ,X 


2 


I  s»(/)dx  =  o.oo?i8  +  J  fi(/)ax  . 

fi  'n  nnnnoA 


0.000026 

The  last  term  in  this  equation  can  be  evaluated  by  a  numerical 
method  because  it  is  a  proper  integral. 

Prom  equation  (33),  we  see  that  U  is  a  function  of  R  for 
any  cross  section  in  the  entrance  region.   If  we  let 
(TJ-3  -  1)R  =  F(R),  then 

J       (u3-l)RdR  =  (      F(R)dR  . 
'0  '0 

This  integral  can  be  evaluated  by  Simpson's  rule.   Thus .the  total 


10$ 


pressure  drop  defined  in  equation  (lj.6)  is  obtained.   Table  2 
(Ref.  1)  shows  the  values  of  the  integrals 


I   0(/)dX  and  f  (U3-l)RdR 


'0  '0 

and  the  results  of  the  pressure  drop  at  different  locations 

from  the  entry. 

Prom  this  table,  we  may  plot  a  curve,  of  against  X. 

1  -2 

2 
This  curve  determines  the  pressure  drop,  p    -p,  as  a  function 

oT  x.   The  slope  of  the  curve  is 

—   ,Po^  ~P}  .      1        ££  _     r02    f£ 

dX  1    r,  1  .  „  dX      1       dx 

-Pu2  -Pu2  (-/tyu) 

2  2  '  2 

In  the  fully  developed  laminar  flow  region  (Hagan-Poiseville 

dp     3P\lu 

hus  the  slope  of  this  curve  at  the  end 


)  =  16. 


2' 
Consequently  the  curve  is  asymptotic  to  a  line  with  slope  16 . 

For  a  value  X-j_,  somewhat  greater  than  the  entrance  length 
Xe,  and  if  X  >   Xe,  equation  ( J+8 }  becomes 


I'll 

JW), 

— 

=   -  

2 

■    .      Thus 

dx 

rC 

of 

the 

er 

itrance 

re 

gion   is 

- 

V 

8/Vu 

{-/>!>  u) 

ro2 
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p_-p        (L     ,         /Al         rK 

^=—  =1  +  2    (u3-l)RdR  +  2     0(/)dX  +  2    0(/)dX  .   (55) 

-f>u2  Xl 

2 

In  this  case,  the  flow  is  fully  developed;  therefore, U  =  2(l-R2), 

and 


A  A  1 

|  (U3-l)RdR  =  8  I   (l-R2)3RdR  -  I 

Yl  In  In 


RdR 


=  8     (R-3R3+3R5-R7) 


dR 


0 


1  1       11 
-  -  3    ■  -  +  3    ■ 

2  4.      6   8 


1  1 

2  2 

Prom  the  definition  of  0(/)  and  the  power  series  for  I  func- 
tions, we  have 

*2  f     1%2 


I2 


PU)   -  —  [ix*  -  i0i2)   =i-  <-±-  -  iQ) 


I2  I2 


$S(iO 


y 


i.y2  +L/k+  . 

!].  16 


1     _      1      ,  L  i     ?      l 

8  96  8  96 


Thus 


/im        (ii  f/  i    _      /im 
y — =>  of"1    "  / — =>  0 


y 


-y2  +L/^+  ... 

8  96 
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i 


■V  + 


Am 


-yk 


-r2  +-A 

8      96 


iim 
/ — y   0  Iq 


J 


=  8  [2  -  l]  =  8  . 
Since  x- — >  &z>   as  / >  0,    It  follows  that  in  the  fully  de- 
veloped region  0(/)  =  8  =  constant.   Therefore. 
,X             ,X 

fit,-/) ax  =  8         ax  =  8(x  -  xx)  . 
xi 

By  setting  I    '/>  (/)  =  K1,  equation  (55)  gives 
'O 

p    -p 

-22 =  1  +  x   +  2K,  +  l6(X  -  X, ) 


. />  uc 


=  16X  +  2K-.+2  -  16X, 


Let 


then 


2Kx+2  -  16X-L  =  m 


P~,  "P 


=  16X  + 


-Pu2 


We  see  from  Table  2  that  if  X]_  =  0.30i+  (/=  O.I4.),  K-,  =  2.5739. 
Thus 

m  =  2  x  2.5739  +  2  -  l6x  0.301).  =  2.28  . 
The  total  pressure  drop  becomes 


10  8 


T3  -P 

00  =   16X   +    2.2t 

1      -? 


61+  x 

+    2.2f 

DRen 
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APPENDIX  III.   SCBLICHTIMG'S  MATCHING  KETHOD 

Consider  that  the  flow  in  a  straight  channel  comprised  of 
two  flat  parallel  plates  has  uniform  velocity  distribution  at 
the  entry.   Furthermore  assume  that  the  flow  is  incompressible 
and  two-dimensional.   In  this  method  (1)  the  flow  field  is 
divided  into  two  sections,  the  upstream  section  and  the  down- 
stream section.   In  the  upstream  section,  the  boundary  layers 
develop  in  the  same  way  as  they  would  along  a  flat  plate  at  zero 
incidence.   The  integration  is  then  performed  in  the  downstream 
direction  so  that  the  boundary  layers  are  calculated  by  finding 
the  acceleration  of  the  central  stream.   In  the  downstream  sec- 
tion, the  velocity  distribution  is  assumed  to  be  the  sum  of  the 
fully  developed  profile  and  a  deviation  velocity  which  approaches 
zero  as  x ^  ©o  .   The  integration  is  then  performed  in  the  up- 
stream direction  in  order  to  find  this  deviation  velocity.   Hav- 
ing obtained  both  solutions  in  the  form  of  series  expansions, 
they  are  joined  at  a  point  where  both  solutions  are  valid.   In 
this  way  an  approximate  description  of  the  flow  field  in  the 
entire  entrance  region  is  obtained. 

Upstream  Solution 

Here  the  x-axis  is  fixed  along  one  of  the  walls  of  the 
channel  in  the  flow  direction,  and  the  ordinate  is  denoted  by 
y1  as  shown  in  Pig.  8.   In  the  downstream  solution,  however, 
the  x-axis  is  taken  as  the  center  line  of  the  channel  and  the 
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ordinate  y  is  measured  from  the  canter  line.  Following  the 
usual  assumptions,  the  boundary  layer  equation  becomes  (see 
Chapter  1,  equation  (11)) 


3u     ^u 

U  —  +  v  — 

5*    dl 


P   dx 


»K 
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(1) 


Let 


£  = 


V>x 


y 


and  assume  the  center  line  velocity,  uQ,  to  be  in  the  follow- 
ing form 

uQ(x)  =  u(l  +  k-]£  +  k2€2  +  .  .  .  )  (2) 

which  is  in  a  power  series  of  V x.   The  coefficients  k-,,  k2, 
...  are  unknown.  '  . 

According  to  Bernoulli's  equation,  the  pressure  gradient 
in  equation  (1)  can  be  written  as 
i  dp      du0 


P   d> 


Mr 


dx 


=  u(l  +  k-j_€  +  k2£2  +  .  .  . ) 


2 

iK^  +  2tk2  +  3izk 


where 


dx    2  '  a2  ux   2  ^  a2  u. 


3  T 

1  t 

2  x 


dx 


Therefore,  we  have 
1  dp   u2  £ 
P   dx    2x 


(1+  k]_£  +  k2£2  +  .  .  .  )  (k 


!  +  2ck2  +  3^23  +  ...) 
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u2e 


2x 


[kx   +    (K12+2K2)£    +   3(K1K2+K3)£.2 


.  (3) 


Introducing  a  new  variable 

n  =  7  —  y 

V  x 
and' assuming  that  the  stream  function,  f[x,    7'),    can  be  ex- 
panded in  a  series  in  £  in  a  form  similar  to  Uq(x)  as  shown 
below 

fix,   y')  =ua[£f0(n)  +  ^(n)  +  .  .  .] 
we  have 


3  f 


u(x,  y')  =  =  ua  j  £f0'(n)  +  £  V(|D  +  •  •  •   — L 

37'  L  d  27' 


where 


Thus 


dr\ 


-i--i 


3  7'  V  x 


1 
a  i- 


u(x,  y')  =  u  [f0'(r|)  +6-f1'()l) 
Similarly, 


(W 


5  r 

v(x,  y')  =  

3  x 


>c24 


-ids 
dx 


-  ua  |fQ  +  2tf1   +   3L*f2   +    . 

-  ua  |£f0'  +  L2tx'    +  L2f2'    +  .  .  .  |  — 


?  x 


Since 


-  u       d1         /u              X       t/o              1 
n   =-{/—  y'  ,   — -  =f  -  y'(-  -  x-V2)  = „ 

Vx       3  x  -J  2  2x 


113 


and 

&l        i 
dx        2x 
velocity  v(x,    y')    becomes 
ua 
2x 


v(x,y')    = I  '  jf  0   +   2ifx  +   3£-2f2   +    •    •  '•]   + 


ua£  n 


2x 


[f0>    +   Cfj.'    +    t2f2'    +    .     .     .] 

Prom  equation    (I4.),    we   have 

^u  r  o  -,   d£ 

—  =  u    if'    +   2tf    '    +   3  erf,'    +    .     .     .      — 
<?x  L1  '*  J  Jdx 


(5) 


|f0"  +  £'1"  +  «-2f2"  +  £3f3"  +  •    •    •  J   — 


tu 


2x 


[fl* 


+   2£fp'    +  3£' 


J3 


]  "  Trl  [V  +  "l" 


42f2"     +    d3f3"     + 


] 


?y' 


=  5  [V  +  £V 


+  t2f0"  +  ... 


11 
?y' 


(6) 


=  —  ff0"  +  ifx"   +   £2f2"   +    . 


d*u       u 


g  7 ' 2      £  a 


f0"'   +  £i\"'+  t2f2"'   + 


1  iT_ 
'  ^y1 


a2£2 


[f0,M    +   £fx'"    +   £2f2"'    +    .    .    .1 


(7) 


(8) 


Substituting   equations    (3)    through    (8)    into   equation    (1),    we 
have 


£u2 
2x 


[f0  '    +    lf± '    +   e2f  2  '    +    •  •  •  ]    If  1 '    +   2£f  2  '    +  3£.2fj  '    +    ...  "I 
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2x 


y\[t0>  +  Lfx'  +  i2f2-  +  ...]§0"  +  ttx"  +  t2f2" 


+  i3f3"  +  .  .  .1 


^   r 


a2£2 


|f0"   +  Efx"   +  £2f2"   +    ...    1  [fQ  +   2£fx  +  3t2f2   +   •••] 

•   ][V    T   «l'    +   i2f2     +    ■•■] 

~2  ■ 
[f0'"    +   Efi"    +  £2f2'"    +    ...]+  ("k-l   +    (K12+2K2)£ 

+   3(K1K2    +  K3)£-2   +    .  ..    I 


2x 

u2 

—  n  [f0"  +  tfy  +  t2f2"  + 

2x 
Vu      r  „  n        u2£ 


where 

1>U  u2 


a2£2    x 
Simplifying,  we  obtain 

t[f0'  +  eiy  +  e2f2T  +  ...  ][fx'  +  2£f2'  +  3£2f3'  +  ...J 

-    |_f0"   +  tfi"   +  £2f2"   +    ...    ]|f0   +  2£il  +  3£2f2  +    •••] 

=   2    [f0'»     +    £fx"'     +  t2f2"<     +    ...]    +f[Kx  +    (Kx2   +   2K2)£. 

+  3(KXK2   +  K3)t2   +   ...1 

Collecting  terras  in  powers  of  i    gives  a  series  of  differential 
equations  for  fQ,  £,,  f2,  ...   .   The  equation  of  the  first 
approximation  for  the  zero  power  of  i.   is 

fof0"  +  2f0"'  =  0  (9) 

with  boundary  conditions  f 0  =  f 0  '  =  0  at  r\   =  0  and  f 0  '  =  1  at 
y\   — >  aO   as  obtained  from  equations  (1|)  and  (5). 

This  is  identical  to  the  differential  equation  of  flow 
over  a  flat  place  given  and  solved  by  Blisius  (1).   The  values 
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L0'  x0  '  10"  can  be  found  in  the  table  (Ref.  1). 

The  second  approximation  (for  the  first  power  of  fc)  is 
determined  by  the  equation 

2lV"  +  f^"  -  f0<tx<    +  2^"^  =  -KX  (10) 

which  is  a  linear  differential  equation  of  third  order.   The 
boundary  conditions  are  £j_   =  f-,  '  =  0  at  f[   =  0  and  f -,  '  =  K-, 
at  y\   =  oo  .   The  first  two  boundary  conditions  are  obtained  from 
equations  {l\.)    and  (5);  the  last  one  is  obtained  from  the  right- 
hand  sides  of  equations  (2)  and  (4),  since  u(x,  /')  =  uQ(x) 

as  Yj >  cx>  . 

For  the  higher  approximation,  we  can  obtain  equations  of 
. irailar  type.   From  these  equations,  values  of  f  and  K  can  be 
found.   Thus^ the  velocity  distribution  in  this  region  is  com- 
pletely defined.   However,  a  great  number  of  terms  would  be 
needed  in  the  series  given  by  (1|)  if  the  entrance  region  problem 
is  to  be  solved  merely  by  this  upstream  solution.   Each  added 
term  would  make  the  solution  more  and  more  difficult  and  complex. 
Therefore, another  mathematical  analysis  is  applied  in  the  down- 
stream region  where  the  deviation  velocity  from  fully-developed 
flow  is  found. 

Downstream  Solution 

Assume  that 

u(x,  y)  =  uf(y)  -  u-*(x,  y)  (11) 

where 


3 

y2 

-  u(l  • 

■  —J 

2 

a2 
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uf(y)  =  -  u(l  -  — )  (12) 

2       a2 

the  fully-developed  velocity  distribution,  and  u-*(x,  y)  is  a 

deviation  of  velocity  from  the  fully  developed  one.   Here,  y 

is  measured  from  the  center  line  of  the  channel. 

Substituting  equation  (11)  into  (1),  we  have 

<ru#      duf   i?u-"- 

(uf  -  u-::-)  (-  )  +  v( ) 

2x       dy    3  7 

d2u7-   ^2u«-    1  dp 

■  V( ) .  (13) 

dy2     ,?y2    P   dx 

Differentiating  equation  (13)  with  respect  to  y,  the  pressure 

term  is  eliminated  because  it  is  a  function  of  x  alone,  as 

follows. 

duf  ^u-»-      ^2u---   ^u»  ju*      j2u-::-  3y     duf   <m* 


v. 


t 


dy  3x  <?x?y   ?y  5x       ?x^y   ,:  y  dy    <?y 

d2Uf   ,?2u*      dJUf   ,?3u» 
+  v( )  =^(— : r-J   •  (1W 

3j2  3J2  31i  31i 

Prom  equation  (11),  we  have 

du  iu-«- 

—  = .  (15) 

3x  ^x 

?u        ^v 

Using   the   continuity   equation,   —  +  —  =0,    equation    (15) 

<?x       31 
becomes 

S  v       0  u* 

—  =  .  (16) 

3  J      \?x 

Substituting  equation  (16)  into  equation  (ll;)  and  simplifying, 

we  obtain 


Uf  +  u-«- 


v  (- 


d2Uf   ?2u« 


?wy 


5x3y 


dy2    py2 
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d^Uf.   ,?3U-;:- 
)  =  v  (— ± -).  (17) 


ay 


?yJ 


For  the  deviation  velocity  u-»-(x,  y)  ,  Schlichting  assumed  a  series 
of  the  form 


n* 


=  C^e  1        ^i  (-)  +  C2e 


r2'(-> 


.y.  -x?e£ 

+  C2e 
u  a 

In  this  series,  only  the  first  term  is  used;  the  higher  order 

terms  are  neglected  in  the  first  approximation.   Thus 


u  a 

Prom  equation  (18),  we  have 


u  e 


Vi"(-) 


3j  a 

5x 


=  u  c1  y1'(-) 


y 

a 

y 


e'Xli    (-2Xl£) 


xie' 


d£ 
dx 
£ 


"J  C-,  tV(-)  e  "-1"   (-2X-,£)  — 


2x 


uXlCl£2   -Xl<2 


fc     rn'(-) 


2„* 


e^u 


jxjy 


uA-iC 


p  2 


ax 


•2         y 

/  » 

l 


e"Xl£'   «S"R 


-    x      rf'(-) 


^y£ 


2y- 


»2 


fxlv(-) 


(18) 


(19) 

(20) 
(21) 
(22>> 
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Prom  equation  (12),  we  have 


aUf. 

ay 

d2uf 

ay^ 


ay 


3 


3 

I2 
3 


o. 


uy 


Prom  equations  (16)  and  (19),  we  have 
su# 


V  = 


ay 


uX1C1e' 


2fy-(-)dy 

a 


.7 


Since 


f-  =M 


3uAl°l£     -Ax£2 

e  x     r  (-) 

x  a 


l»x 


(23) 
(22j.) 


(25) 


\)  = 


£2a2u 


(26) 


Substituting  equations  (12)  and  (id)  and  equations  (20)  through 
'26)  into  equation  (17)  yields 


3 


uX1C1£' 


-  u(l  -  — ) 


e-*l£  n" 


uA^f 


*V 
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su\-,C-,S^        ,    ,2  3  uC-,        ,,£2 

iJ_  s-^l£     r,( u  -  — i  e-^     r,"') 


Id  ? 

a'  n 


uCl       .\j£Z 


.2  =  2-. 

(0 


x 


Simplifying   the   above   equation,    we   obtain 

u  2  a^  J 

=  C1\1e"Xie      (  4V   n"   "  ^i  Vi"')    • 

The  right-hand  side  of  the  above  equation  is  zero  as  higher 
order  terms  of  •/-,  can  be  neglected.   Then  we  obtain  the  follow- 
ing equation. 

VniV  +  3*1  |¥,  +  -  (1.  -  ^r)  Vl  =  0  (27) 

with  \-,  as  it3  eigenvalue.   The  boundary  conditions 
v  =  0  and  —  =  0  at  y  =  0 

■ay 

v  =  0  and  u  =  0  at  y  =  a 

are  equivalent  to 

:,x  =  yx"=   o   at  y  =  0  (28) 

vl   =  ^i'  =  0  at  y  =  +  a  .  (29) 

The  other  boundary  condition  is  chosen  as 

yx'  =  1   at   y  =  0  .  (30) 

This  is  possible  because  C]_  is  still  free.   Having  these 

boundary  conditions,  the  equation  (27)  can  be  solved  by  a  power 

series  method. 
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Let 


VjH  ■=  L    An(-)n 


n=0 


An  +  An (-)  +  Ao(-)   + 


a       a 


Then  we  have 


y,n-i 


V  =  ^  nAn(-) 
n=l 


y        y  2 
=  A1  +  2A2(-)  +  3Ao(-)   +....- 
a        a 

^"  =  r   n(n-l)An(-)n-2 
n=2  a 


y        y  o 
2A2  +  6A3(-)  +  12A^(-)   +  .  .  .  . 


Yi"    =  Z     n(n-l)(n-2)An(-)n"3 
n=3  a 

y         y  p 
=  6A3  +  2i+Ak(-)  +  oOA^(-)   +  .  .  . 
a         a 

^1V  =  T  n(n-I)(n-2)(n-3)An(-)n^ 
n=4.  a 

y  y  o 

=  21|Al  +  120A^(-)  +  360A6(-)   +  .  .  .   . 
a  a 

Substituting  the  above  series  for  y-,  ,  f,  " ,    and  m/1V  into 

equation  (27),  we  obtain 

ZT  n(n-r)  ,n-2)(n-3)An(^)n"^ 
n=l+  a 
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3*1 


<■       y  n     1         y2    ^  y  n-2 

1  An(-)n  +  -  (1  -  — )  2_  n(n-l)An(-)n 
n=0    a     2      a"2   n=2  a 


V  f 


n=0 


(n+^)(n+3)(n+2)(n+l)An+lj_{-)n 


3*i 


^      y      1  <*?  y 

r  A(-)n  +  -Z_   (n+2)(n+l)An+2(-)n 
n=0    a     2  n=0  a 


1  r-  y  n 

-  -  £_   n(n-l)An(-) 


2  n=2 


Therefore .we  can  write 


=  0  . 


r     1 
lj.JA^  +  3X1  |_A0  +  -  (2-l)A2  J  =  0  , 


5.^.3. 2A5  +  3\1  ^A1  +  -  (3-2)A3 J  =  0y 


and 


(n+i|.)  (n+3)(n+2)  (n+l)An+1, 


+  3^n 


1-1 

-  n(n-l)jAn  +  -  (n+2) (n+l)An+2  \     =  0 


The  general  recursion  formula  for  the  coefficient  then  is 
_  3^x  [(n-2)(n+l)An  -  (n+1) (n+2) An+2  ] 


(3D 


2(n+l) (n+2) (n+3) (n+k) 
From  the  boundary  conditions  shown  in  equation  (28),  we  obtain 
Aq  =  A2  =0.   Therefore  Al  and  all  other  constants  of  even  sub- 
scripts are  zero  as  seen  from  recursion  equation  (31) .   All 
constants  with  odd  subscripts  from  five  to  infinity  can  be 
represented  in  terms  of  A]_  and  A3,  therefore.  A-^   and  A,  are  left 
as  arbitrary  constants.   Prom  the  boundary  condition  given  in 
equation  (30),  A]_  is  found  to  be  1.   The  remaining  two  boundary 
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conditions  V1(+  1)  =  0  and  f]_'(±  1)  =  0  as  shown  in  equation 

(29)  are  used  for-  determining  Aj  and  the  eigenvalue  7u.   Since 

these  boundary  conditions  are  in  the  form  of  an  infinite  series, 

A-j  and  \-^   are  found  by  a  trial  and  error  method.   First,  a  value 

of  X-j^  is  chosen,  then  different  values  of  An  are  found  in  terns 

of  A^  by  equation  (31).   Using  one  of  the  remaining  boundary 

conditions,  A,  is  found.   Then  the  remaining  boundary  condition 

is  used  to  check  the  values  of  A^  and  X-^.      If  these  values  do 

not  satisfy  the  last  boundary  condition,  a  new  value  of  X]_  is 

chosen  and  the  process  is  repeated  until  the  correct  value  of 

X]_  is  found.   Schlichting  gave  the  value  as  \\   =  18.75.   In  this 

manner,  all  the  constant  terms  can  be  found.   Substituting  these 

y       y 
constants  into  the  assumed  series  forlf1(-),  'A'(-)  is  defined 

a       a 

y 

■  various  values  of  —  . 
a 
For  the  complete  solution  of  tne  velocity  profiles,  only 

Cm  and  the  values  of  £   at  the  joining  point  remain  to  be  found. 

These  are  founa  by  matching  the  upstream  and  downstream  solutions. 

Upstream  Press-ire  D5  -tribution 

Near  the  entrance  of  the  channel,  the  boundary  layers  are 
very  thin.   Most  of  the  flow  is  composed  of  the  accelerating 
potential  flow  in  the  central  core  of  the  channel.   Therefore, 
the  pressure  distribution  is 
j.  dp      dug 
?   dx      dx 
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Integrating  from  x  =  Xq,  where  p  =  pQ  and  uQ  =  u,  to  any  section 
in  the  upstream  region  gives 


i  /?     1  fuo 
1   PO         u 


P  -  P0  =  -  />(U2  -  uQ2) 


P  -  P0    ,  ,u0,2 

=  1  -  (  — ;   .  (32) 

-  f>   u2 
2 


Downstream  Pressure  Distribution 


From  equation  (13),  we  have 

1  6p    ,,d2uf   ?2u*              3u-::-     duf   ju-:<- 
-  -  =  V(  —  -  — )  +  (Uf  -  u») v{-i ) 

Y   ex      ay<^     <Jy2  2X  dy     ^y 

1  dp 
Since  the  pressure  varies  with  x  only,  can  be  evaluated  at 

P   dx 
y  -  +  a.   Therefore  all  terms  except  the  first  two  on  the  right- 
hand  side  of  the  above  equation  vanish.   The  equation  then 
becomes 

■1  dp     d2uf   ^2u-tf 

~p7-  =  iJ{—2 T)  ^3) 


Substituting  equations  (21)  and  (23)  into  equation  (33), 


we  nave 


1  dp     f      3  uC-,       2        -r 
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It  is  to  be  noted  that  Y-^"  '  (-)    Is  defined  for  all  values  of  -  , 

a  a 

J  /j>x 

since  V-,  (— )  has  been  found.   By  definition  we  have  £  =y — — 

a  a2u 

Therefore 

dx  =  2aRea£.d£  . 
Let  the  vsxue  of  £  at  joining  point  be  denoted  by  £-:;-.   Inte- 
grating equation  (34)  from  £  =  £.#,  where  p  =  p-::-,  to  £  ,  one 

.obtains 

P       3vu         f1      DuC-,  fc'      2 

dp  = •  2aRea   £d£ =  w  '"  (1)  -2aRea   e"Al  d£ 

p#  a2  'e*  a2  '£* 

p  -  p*      6ReQ    „    „     2C-,   Re.  ,nir2      »,f,2 

1  „      ua  X-i   ua 

2  y  " 

or 

P  -  P*        o    o    2C-,  .  e2  ,  ,,,2 

=  .  6(t2-  £*2)  +  — i  f  "'(l)(.e"Xle   -  e"Al£"  )  .   (35) 

2  ! 

For  the  complete  solution  of  the  pressure  distribution,  only  C, 
and  £.*  remain  to  be  evaluated.  These  values  are  determined  in 
the  process  of  matching  the  upstream  and  downstream  solutions. 


Matching  of  Upstream  and  Downstream  Solutions 

The  upstream  and  downstream  solutions  of  velocity  and  pres- 
sure distribution  are  to  be  matched  at  a  point  where  both  solu- 
tions are  valid.   Because  there  are  two  undetermined  quantities, 


Uf.(0)  -  uC-g"Xl£'   V,  '(o) 
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C-^  and  £-"-,  in  the  equations  to  be  joined,  two  conditions  are 
needed  to  solve  for  them.   One  condition  is  that  the  center  line 
velocities  of  both  solutions  are  equal  at  the  joining  point. 

Another  condition  used  is  that  the  slopes  of  pressure  distribu- 

dp 
tion  —  ,  for  the  two  solutions  match  at  the  joining  point. 

dx 
These  two  conditions  are  used  to  develop  two  simultaneous  equa- 
tions from  which  Cx  and  £""'  are  determined. 

The  downstream  center  line  velocity,  from  equations  (11) 
and  (18),  is 

Uq  =  u^ioj  -  uU- e  "~J-   y,  ' 
Because  of  equation  (12)  and  the  boundary  condition,  equation 
(30),  the  above  equation  can  be  written  as 

3_   _    -x-t? 
u0  =  —  u  -  uC-j_e   J-    .  (36) 

Equating  the  center  line  velocity  represented  by  equations  (2) 

^nd  (36)  gives 

3  2 

1  +  Kxe   +  K2iA    +  K3£3  +  \ihr  = C1e"A'l£   .  (37) 

It  is  to  be  noted  tha   the  infinite  series  of  equation  (2)  is 
terminated  at  the  fifth  term  because  only  K, ,  K?,  K,,  and  Ki 
.iave  been  four./.  (Kef.  2). 

The  upstream  pressure  slope  equation  is 

dp         dug 

—  =  "  f*o   —   • 

dx         dx 

Substituting  the  center  line  velocity  in  the  upstream  section 
gives 
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—  =  -  fv.2{l   +  Kx£  +  K2£2  +  K3£.3  +  Kl£^)(Ki  +  2K2t  +  3K3£2 
dx 

,  d£ 

+  1|Kl£j)  — 

-    (1  +  Kx£    +  K2^2   +  K3e3   +  K^)  (K^   +   2K2i    +   3K3^2 


2iaRe 


a 


+  ^3)  .       (38) 
Equating  equation  (31+)  ,  the  pressure  slope  equation  in  the  down- 
stream section,  and  equation  (33),  we  have 
3»u   VUG-,  .\-,c2 


,2 


e    -     n' 


(1  +  Kx£    +  K2£2   +  K3£3   +  K^)  (Kx  +   2K2t 


or 

:2  1 


2£a  Rea 

+   3K3^2    +   1+K^3) 

■\-,L2  1 

3   +  Cie"A1        V1'"(l)    =  —   (1   +  Kx£    +  K/2   +  K3&3   +  K^t^-) 

(Kx   +   2K2t  +   3K3£2   +  1+K^3)  (39) 

Eliminating  C^e  from   equations    (37)    and    (39)    gives   a 

seventh  order  equation   in   L    as   follows    (Ref .    2) . 

Kl  +  ^[Kl2   +   2K2   -   6    -  V1"'(l)]    +   ^2[3K3_K2   +   3K3   +   ZZtfjt'llj] 
+   ^[l+K-jKj    +    2K22    +  1+K^  +    2K2yi"'  (1)]  +  £^[5KXK^  +   5K2K3 

.      +   2K3V1"'(1)]    +   £5[6K2K^  +  3K32   +   2Ki|^1"'(l)] 

+   TKjK^6   +  IfK^V     =  0  (40) 

The   values  of  K ' s  were  not  given  in  Sohlichting ' s  paper.   In 
Ref.  (3)  they  are  Kj_  =  1.7207,  K2  =  -2.233,  K3  =  19.81,  and 
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Kr  =  -208.896  (2).   Schlichting  gave  the  value  I   at  the  match- 
ing point,  that  is,  £  =  t*,    as  0 .  126I4.,  and  substituting  this 
value  into  equation  (37)  yields  C^  =  O.3I4.85. 
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For  many  years,  the  study  of  the  developing  velocity  pro- 
files and  pressure  drops  in  the  hydrodynamics  entrance  region  of 
tubes  and  ducts  has  attracted  the  interests  of  many  investi- 
gators.  Although  the  exact  analytical  solution  is  not  available 
at  the  present  time,  many  approximate  solutions  have  been  pre- 
sented.  The  approximation  methods  of  solution  can  be  classified 
into  four  different  categories:   the  momentum  integral  method, 
linearizing  method,  matching  method,  and  finite  difference 
method. 

The  purpose  of  this  report  is  to  present  a  detailed  study 
of  these  four  different  approximation  methods  and  to  present  the 
results  obtained  by  applying  a  finite  difference  method  in  solv- 
ing an  MED  entrance  region  problem  in  which  the  velocity  profiles 
are  developed  from  a  nonmegnetically  fully  developed  one  (para- 
bolic) at  the  entry  to  a  Hartmann  velocity  profile. 

A  brief  introduction  of  the  Navier-Stokes  equation  and  the 
boundary  layer  equations  is  given  in  Chapter  1.   In  Chapter  2 
Schiller's  solution  by  the  momentum  integral  method  for  a  flow 
in  a  circular  tube  is  presented  in  detail.   A  momentum  integral 
equation  is  derived  first.   Then  by  applying  the  assumed  velocity 
profile,  which  involves  the  boundary  layer  thickness,  to  the 


mo: 


nentum  integral  equation,  a  relation  between  the  boundary 


layer  thickness  and  the  axial  coordinate  is  obtained.   The  solu- 
tion is  completed  after  carrying  out  the  integration. 

Cnapter  3  contains  a  detailed  account  of  Langhaar's  lin- 
earizing method  for  the  flow  in  a  circular  tube.   The  nonlinear 
partial  differential  equation  of  motion  for  the  axial  direction 


is  linearized  by  the  assumption  that  the  conveotive  terms  are  a 
function  of  axial  direction  only.  The  solution  is  presented  in 
closed  form  in  terms  of  Bessel  functions. 

Chapter  1+   concerns  Schlichting' s  matching  method  in  solving 
the  entrance  region  problem  of  a  flow  in  a  parallel  plate 
channel.   The  flow  is  divided  into  two  sections,  the  upstream 
and  downstream  sections.   The  solution  for  the  upstream  section 
corresponds  to  the  Blasius  series  solution  over  a  flat  plate, 
and  the  solution  representing  the  progressive  deviation  of  the 
profile  from  its  asymptotic  parabolic  distribution  is  obtained 
in  a  series  form  for  the  downstream  section.   Solutions  for  these 
two  separate  sections  are  then  matched  at  a  point  where  both 
solutions  are  applicable. 

In  Chapter  5,  Hwang  and  Pan's  finite  difference  method  is 
applied  in  solving  the  MHD  entrance  region  problem  with  the  new 
boundary  conditions  mentioned  above.   The  continuity  equation 
and  boundary  layer  equation  are  reduced  to  finite  difference 
equatio:. s  and  solved  on  a  digital  computer.   The  results  ob- 
tained are  compared  with  those  obtained  by  a  momentum  integral 
method. 


